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Abstract 


In this dissertation the problem of PN code acquisition foi DS CDMA. sa steins 
operating m low SNR environments has been investigated Among the vinous tedi 
mques Stiial seaich code acquisition technique has been us< d md impim eel upon The 
modified stiategy has been analyzed from fixed and adaptive threshold detection point 
of view in unfaded incl faded environments chai actei ized In Nakiguni m f ichng Im 
poit nil acquisition parameters namely mean ( T c q ) and \ unnee (<r^ ( „) of acquisition 
tune have been derived using the signal flow graph technique The results obt un<cl 
demonstrate that foi low SNR improvement of one oidei ol magnitude is possible m 
the aloiesaid acquisition parameters when compared with existing techniques 
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Chapter 1 


Introduction 


1 1 Introduction 


A. Spread Spectrum (SS) system produces a signal whose emigy is distributed oui i 
bandwidth much wider than the message bandwidth As such it is characterized bv l 
bandwidth expansion factor N defined as the latio of bandwidth required to transmit 
a SS signal (W) to the message signal bandwidth (B) i e 




( 1 1 ) 


The most commonly discussed spread spectrum systems are Direct Sequence (DS) 
and Frequency Hopping (FH) [2 4, 6 9] Less commonly discussed are chnp md 
time hopping [9] Additionally these techniques can be combined to present a hybiicl 
system which combines the attractive features of individual techniques Of the numbci 
of possible techniques, we restrict our attention to Direct Sequence Spread Spectrum 
(DS SS) 

Figure 1 1(a) is a simplified block diagram of a DS SS transmitter [3] The bise 
band data waveform is first converted into the binary signal Ihe resulting symbol i itc 
is denoted as A- The next operation is multiplication with the binary pseudorandom 
(PN) sequence whose chip rate is denoted as jr The PN code chip rate is chosen 
much faster than the symbol rate so that T s = NT C As illustrated m Fig 1 1(a) (he 
effect of this operation is to “ spread ” the bandwidth of the waveform by a factoi N 
which for the same signal power causes the spectral density of waveform to be quite 
low and noise like The signal power as indicated in Fig 1 1 ( i) is reduced by a factor 
of A from its level without spreading 

At the receiver, “ despreading ” (multiplication by the same binary antipodal spread 


1 




Received 

Spectrum 


(b) Receivei 


Figure 11 (a) A simple 

Receiver [3] 


rum transmitter (b) 
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mg waveform as at the transmitter) is performed first and then removal of the carnei 
modulation is done to restore the original base band data waveform allowing the le 
ceiver to filtei out a large part of wide band interference Assuming as m Fig 1 1(b) 
the receiver front end filter admits the W hertz wide desired signal The interference oi 
jamming within this bandwidth (BW) will also be admitted b\ the receiver Howeiei 
level of interference (denoted as A/o m Fig 1 1(b)) can be consult rably laiger than the re 
ceived signal (denoted as 4 X ), implying an RF signal to noise i \tio (SNR) rf = At < 1 
But after despreading the bandwidth of desired signal is icduced to its oiignnl a il 
ue B while the bandwidth of the interference remains W Thus filteimg to the signal 
bandwidth can be used to eliminate interference power resulting m the base band data 
SNR 

(SNR) basebond = = Njjr = ^ (SNR) Rr ( 1 2 ) 

Therefore DS SS is a method of using a periodic PN code to generate a large B"V\ 
signal having considerably low power spectral density (PSD) as compared to conven 
tional systems Unlike most other signalling schemes or communication systems which 
are typically bandwidth limited, spread spectrum is conveisc 

Historically DS SS has its roots in the military arena while desirable commumca 
tions lequire Jammer resistance, Low Probability of Intel ccpt (LPI) and Low Proba 
bility of Detection (LPD) and Low Probability of Exploitation (LPE) as discussed m 
[7, 8, 9] However DS SS has gamed interest in civilian communications commumtv as 
the basis of the multiple access scheme known as Code Division Multiple Access (CD 
MA) A well known commercial system is IS- 9 5 (ongmally designed by Qualcomm) 
Many chip sets aie being developed m support of CDMA not only by Qualcomm but 
also Stanford Technologies and Loral (paramox) [35] 


1 2 Synchronization in Spread Spectrum System 

In a typical spiend spectrum system, any meaningful transmitted information can be 
extracted from the leceived signal if 

• Code at the leceiver is accurately aligned with the spi( iding code contained m 
the received signal i e, the two codes must have identic il chips and then edges 
(timing) must piecisely coincide 

• The transmitted and receiver carrier are accurately aligned 
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If the two codes are out of alignment by as little as one chip insufficient signal eneig\ 
will reach the receiver data demodulator for reliable data recovery The process of 
piecisely aligning the two codes is termed as Code Synchronisation [1 2] It pins i 
cardinal role in the proper functioning of any spread spectrum system The piotc ss of 
code synchronization is performed in two steps namely code acquisition and code 
tracking Code acqmsition is a process of successive decisions wherein the ultnn ite 
goal is to bring the two codes into coarse time alignment within one code chip intc n il 
After this condition is detected and verified, the tracking systc m is activated Tuicl mi j 
also called fine synchronization is the operation by which suiclnomzation enois aie 
further leduced or at least maintained within certain bounds [9, 20] Coupled with 
the code synchronization issue is the synchronization of transmitted carrier with that 
of receiver carrier The code and carrier synchronization aie inter linked and one can 
not be achieved without the help of other Therefore, following steps aie imohed foi 
achieving the two synchiomzations 

1 Coaise carrier frequency estimation to get an idea of nominal carnei frequency 
in use 

2 PN code acquisition 

3 PN code tracking 

4 Coarse carrier phase estimation 

5 Fine alignment of carrier fiequency and phase 

Thus, PN sequence timing can be acquired and tracked accur ately and hence spreading 
can be lemoved without the knowledge of carrier phase and with only rough estnn ate 
of earner frequency Once code timings has been acquired, earner phase and fiequuua 
can be accurately estimated using the conventional phase locked loop techniques in 
the same way as for any digital communication system [1] 

This thesis deals with the first part of code synchronization namely ' Code Aiqut 
sition” 


1 3 The Acquisition problem 

In a DS SS systems, the complexity and random nature of initial code acquisition pio 
cess is due to a number of factors such as 1) initial uncertainta about code phase offsc t 
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2) channel distortion (e g fading channels) and additive inter k lence (intentional oi un 
intentional) 3) unknown carrier phase and possibly carrier fiequency offset (Dopplei) 
4) front end receiver additive white Gaussian noise (AWGN) 

The set of design parameters for the acquisition procechuc include threshold sit 
tings correlation time i e number of test per code chip and system complexity as 
manifested by the choice of search strategy, verification logit etc Also implicit is the 
knowledge of important parameters such as design SNR code l ate code length code 
uncertainty legion etc As a result of the above factors th e time to acquire the code 
also called acquisition time ( T acq ) is a random "variable and its complete statistical 
description requires probability distribution function Although the pioblem of de 
termming the distribution function of acquisition time can be solved in principle it is 
very difficult to obtain in practice at least in closed foim However, as pointed out in 
[2 6, 20] the fiist two moments namely mean acquisition time T acq and variance 
of acquisition time o 2 acq provide sufficiently accurate statistical description for most 
practical pui poses 


1 4 Aim of thesis 

The code acquisition problem in spread spectrum systems Ins attracted consider able 
research in rece rt past [15 16 17, 18, 20] As shown in fig 2 2 a number of viable sc rich 
strategies are possible for code acquisition Prominent among those are maximum 
likelihood (ML) appioach, whose parallel implementation complexity is prohibitive 
Howevei a less complex senal version for ML approach is possible in which a decision 
can be made only after examining the outcome of each cock phase On the otlici 
hand low SNR environments are characterized by stiong lnteifuing signals, deep noise, 
Doppler effect and rapid fading etc therefore multiple examination of each code pli ise 
m serial version of ML approach can not be ruled out [20] Sequential approaches 
are coheient m nature and have lower implementation complexity than ML method 
Then merits strongly depends on the accuracy of the initial estimates of code phase 
which reduces quite rapidly with decreasing SNR and then foie, not well suited nr 
typical low SNR environments [6] The serial search approach on the other hand is 
charactenzed by its simplicity, reliability, strong interference rejection capability m 
hostile cmnonment, ability to accommodate long PN sequences and low hardware 
requirements Because of its attractive features the majority of research work on PN 
code acquisition revolves around the serial search strategy 

Recently [15] has proposed a very attractive serial seaie h strategy foi low SNR 
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CDMA systems based on the recirculation loop punciple m which the main reason hn 
low SNR was cited as considerable code Doppler and receix c d signal was assumed to 
be unfaded The analytical frame work of [15] was pure prob rbihstrc which except foi 
the cases can not be generalized The dynamic SNR range c onsidered there w is not 
sufhcently wide from spread spectrum point of view and m that range another method 
namely sequential estimation has better performance m teims of code acquisition pi 
rameteis Moreover to improve the detection of correct code phase the constant, f list 
alarm probability approach was utilized using a complex Ad iptive Decision Thushold 
Level Conti ol (ADTLC) [17] As such the work of [15] has left many open issues which 
require careful addressing Some of these issues forms the bisis of work embodied m 
this thesis with the aim to modify the scheme of [15] on two fronts 

1 The addition of post detection mtegiation feature to widen the dynamic SNR 
lange so that lelatixely simple hardware could be used foi code acquisition 

2 Modification in decision cncuitry to achieve the const mt false alarm probxbihtx 
by introducing Cell Averaging Constant False Alarm Rate (CA-CFAR) ap 
proach which requires relatively low hardware as complied to ADTLC method 

These modification will introduce several new issues and some of them will be an xlx /eel 
in this work in particular, our attention will be on 

• To provide an analytical frame work for unfaded recen c cl signal for the modified 
senal search strategy fiom detection point of view 

• To peifoim the analysis of modified senal search stiategx foi low SNR opei it urn 
Here the case of low SNR due to rapid fading is chosen For analytical fi tint 
work the Nakagami m fading model is used which is better suited for modeling 
of urban multi path radio channels [32] 

• Use of more powerful and versatile signal flow graph bast cl approach w ill be m idt 
to derive the geneial form of various acquisition paranuteis 

• Finally, the issue of how much improvement the modified scheme can pi ox ldo ox < i 
the known serial search scheme including the one piesented in [15], for impoit uit 
acquisition parameters m faded and unfaded environment will be considered 

To summarize, the main motivation for this work is to find new serial search alqo 
nthms suitable for all digital receiver implementation and operation at low signal to 
noise ratios 
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1 5 Thesis preview 

In chapter 2 a review of code alignment methods is presented In chaptei 3 the mod 
lfications to senal search scheme is earned out and analyzed fiom detection and code 
acquisition point of view The unfaded and Nakagami m faded signal case is consul 
ered separately The detection and false alarm probability expiessions for the two c isc s 
have been derived Finally using the signal flow graph appio icli the impoitant rcqui 
sition paiameters are derived The detailed derivation imohed m each of the lbou 
case is not included m the mam body of chapter 3 rathei they are pio\ided in the 
appendices In chapter 4 some of the numerical lesults ai< pnsented and discussed 
Finally in chapter 5 scope for further investigations alongwith concluding l email s are 
presented 


Chapter 2 


PN Code Synchronization in Direct 
Sequence Spread Spectrum 
Receivers 

2 1 Introduction 

This chapter examines the following areas which influence s-\ nchronization of DS SS 
systems spieadmg codes coarse synchronization and fine sjik hi onization Spieading 
codes are the highei rate binary sequence of symbols which are used to spread the 
data sequence Coarse synchronization or ‘Acquisition and hue synchronization or 
‘Tracking” together complete the synchronization process of uceived and localh gen 
erated PN code Duung coarse synchronization the most important issue is usualh the 
minimization of acquisition time while SNR maximization and hence minimization of 
bit eiroi (BER) is the issue of concern during the fine svncln onization 


2 2 Literature survey 

If one searches through the available PN code acquisition lilaatuie the common de 
nominator found among almost all the methods is that lecemcl and local PN signals 
are first multiplied to produce a measure of correlation betwc ( n the two This coneh 
tion measure is then processed by a suitable detector/decision rule and search stiatcgr 
to decide whether the two codes are in synchronism and wh it to do if they aie not 
The difference between the various schemes depend on 
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• Type of detector and decision strategy used 

• The nature of search algorithm which acts on the detector output to reach on 
the final decision of code phase acquired or not 


All known detector for PN acquisition systems fall into two basic categories nameh 
coher ent and non coherent By far the most common found in acquisition sj steins foi 
DS SS receivers is the non coherent detector which foi example might be compuscd 
of a band pass filter centered at the frequency of received camel upon which the PN 
signal is directly modulated followed by a square law envelope detector an mtegi ite 
and dump circuit and a simple thieshold device [34] 



Figure 2 1 Structure of detectors for code acquisition [20] 

Consideration has also been given m the literature to PN requisition systems which 
employs a coherent detector, typically implemented as an integrate and-dump ennut 
followed by an optimum Bayes detector [6] or instead just i simple threshold dewu 
[25] The inherent assumption here is that receiver is capable of determining good 
estimates of carrier phase and frequency offsets 
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Fuithei classification of detectoi is based on the type of liihgiation time It ma\ lx 
fixed 01 variable Under the category of fixed integration turn another subdivision is 
possible namely smqle dwell time [23 34] and multiple dwell time detectoi [22] Tlx 
operation of former is based on single fixed time observation of received signil plus 
noise while latter considers multiple such observations Depending upon the duration 
of observation or equivalently the time alloted to make i decision rehtive to tin 
PN code penod, single dwell time detectors can be fuithei dilh lcntiated according to 
whethei thej utilize the full penod oi partial period code emulation The multiple 
dwell detectors differ from one another in the way in which tlx additional observations 
aie used to verify the temporary decision made based on tlx fust observation alone 

Finally the categories of variable integration time is reserved for those cases where 
the dwell time (time for continuously integrated stochastic prore ss to exceed a threshold) 
is a landom variable Typically they employ the classical method of sequential detection 
[2 4] 

The next classification of detectoi stiuctuie is in accoidanu with the i ite at which 
the decisions are made on each code phase position undei test High decision late 
detectors, such as used in matched filter (passive conelation) P\ acquisition sj stems 
[6, 21], lefei to those stiuctm.es that make their decisions it a PN code chip i ate oi 
integer multiple of it Low decision rate detectors which eniplov active coirelition 
make decisions at a rate significantly slower than code chip l ite 

To complete the detector classification further categonz ition is based upon tlx 
criterion used for deciding between in sync or out of sync hvpotheses eg Boijps 
Neyman Pearson and Mini max criterion 

A summaiy of above classification of detector stiuctures used for PN code acquisi 
tion is illustrated m Fig 2 1 

The remaining aspect of the acquisition process is the sc iicli sti itegv by which 
one denotes the particular procedure adopted by the receivei in its seaich tlnough the 
uncertainty region Not unlike the classification of detectoi tv pes the classification of 
acquisition scheme accoidmg to seaich strategy follows a tic c stiuctuie [20] as shown 
m Fig 2 2 Conceptually the simplest of search techniques on the first level ol tlx 
hierarchv is the maximum likelihood algorithm It requires tint leceived PN sign'd is 
correlated with all possible code positions of local PN code u plica The coirect P\ 
alignment is chosen as that local code phase position which pioduccs the maximum 
detectoi output 
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Another scheme is based on a sequential estimation of tin shift register states of 
PN generator typical example is RASE (Rapid Acquisition In Sequential Estimation) 
technique first introduced by Ward [25] In this technique the best estimate of fust 
n chips of leceived PN code is loaded in the n stage local PN code generator thus 
defining a particular initial condition from which the generatoi begins its operation The 
success of this scheme depends on the accuracy of initial cstim ite 

The senal synchronization of PN sequence was fust introduced b\ Sage [24] In its 
ongmal version a senal search was performed by linearly \ uwng the tune diffeunct 
between the PN modulation on the incoming carrier and the PN waveform generrted 
at the receiier with a continuous decision process determining when synchronization 
is achieved Such a system is also called sliding correlator [4] Sweeping through the 
time uncertainty region can also be done in disci ete steps win k the time unceitunh 
region has been quantised into a finite number of elements callc cl cells through which 
the receiver is stepped Which particular search stiategy is sc lccted by the receivu is 
dependent on the nit uie of uncertainty region available pnoi information statistic il 
quality of the tests performed, availability of stepping and n winding mechanisms \ 
multitude of such search algorithm is shown in Fig 2 2 If the uncertainty region co\ pis 
the whole code penocl a stiaight line search would b( appiopnate [19 24] On the 
other hand pnor information is sometimes available which suggests the use of moic 
sophisticated strategies ( Z search or Expanding window searih) for bettei results [19] 
Finally it may be noted all the above serial seaich technique s base then decisions on 
the likelihood ratios ind thus, are fundamentally different fiom sequential estimation 
techniques [25, 26] whose merits decreases with the deci easing SNR 


2 3 Spreading codes 

Direct Sequence (DS) systems require a pseudo random noise (sequence) geneiatoi it 
transmitter and receiver The desirable properties for the spic ichng codes of a single 
transmitter ind receivei pair are 

• An impulse like auto correlation function (oi equivalc lit ly a flat power spectnl 
density) so that when spieadmg and despreading cocks aie not aligned a nc u 
zero coirelation lesult which has numerous implications 

1 The larger the ratio of peak to side lobe of auto coirelation function, the 
easier it is to determine the correct spreading cock ph ise position and lit nc e 
synchronize the despreading sequence with the ret c ived signal 
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Figuie 2 2 Classification of acquisition schemes [6] 

2 Piovides laigcr processing gam and hence bettci utilization of alloc ited 
spectrum lc idmg to increased capacity 

3 can be geneiated m synchronism at transmittei md leceivei 


The piecedmg piopeities are also desirable foi multi usei < ommunications such as 
CDMA additional desirable properties include 

• One spieadmg code pei CDMA user l e sufficient nunilx i of spreading code foi 
all users 

• Orthogonality between spreading codes of diffeient useis so that the extr rction 
of desired usei’s signal could be accomplished from multiple access interference 
(MAI) 

Some of the spreading codes are Maximal Length (m) Secpicnces, Gold codes [29] 
Rapid Acquisition codes [5 28] 
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2 3 1 Maximal length (m) sequences 

Traditionally the maximal length (m) sequences aie attractive because they can bt 
geneiatecl via a Linear Feedback Shift Registei (LFSR) Tilt length and hence ilso 
repetition period of an m sequence constructed fiom an n st ige shift registei is 


L=2 n -1 


(2 1 ) 


Maximal length sequences have a good balance of 1 s and 0 s and when mapped to 
an antipodal signalling scheme have a two valued normalized periodic auto correl it ion 
function 9p(k) given by 


e p (k) = 



(2 2 ) 


wheie l is any integer and L is the sequence period Thu c fore over the internal 


r |< 


LT, 


the auto correlation of maximal length sequence is given as 


h- 

(L± 1 M' 

V L T J 

) 1 T 1 

1 . . 

1 L 


1 T 1 


R p (r) 

and the powei spectial density (PSD) is given as [9] 

L + 1 ^ o /1\ . ( . i 


S P (f) = 


L 2 


E 


sine 


iM' 


LT 


+ T Af) 


l = -CO 
I ^ 0 


(2 3) 


(2 4 ) 


The two functions are plotted m Fig 2 3(a) and (b) The PSD function consists of dolt i 
functions sepaiated by -pp 

As pointed out earlier the more impulse like the auto correlation function of i 
sequence the more preferable it is However, the undesirable pioperty of m sequence 
foi DS SS application is cross correlation between the m seqvc nces can be high affecting 
multiuser separation and code synchronization Another undesirable property of rn 
sequence is that they are inherently susceptible to mathematic al cryptanalj sis Thus 
if message security is desired when m sequence is used for spn admg, the data symbols 
must be enciphered before they arc combined with pseudo noise (PN) 1 sequence [9] 

The poor cross correlation deficiency of m sequence is ic moved m another fannh 
of codes called Gold codes A Gold code is constructed by the modulo 2 addition 
of specific cyclic permutations of a preferred pan of m sequences The period of un 
code in the set of gold codes is L, where L is given by (2 1) The improvement m 


1 In literature the term m sequence and PN sequence are used inter dnngeably 


Chapter 2 


14 



Figure 2 3 (a) Auto correlation function of a m sequence (b) Power spectral 
density of m sequence 

cross correlation between the gold codes in a set, simple const met on of the geneialois 
and abundance makes gold codes attractive foi multiusei communications The auto 
correlation however is no longei as attractive as that of m sequences and consequuith 
the acquisition process can suffer Gold codes and then propeities are considered in 
moie detail in [29] 

2 3 2 Rapid acquisition sequences 

The problem of rapid code acquisition has received consider ible attention Stifflei [28] 
investigated the possibility of designing codes which have lapicl acquisition proper tus 
The principle of operation was that a correlation is performed uid depending upon the 
correlation value, if another correlation is required the uncert unty region of the code 
is halved Thus the average number of correlation to achieve the acquisition is reduced 
to log 2 L from ~, although the correlation time for each step of acquisition pioass is 
longer than that required in a conventional system 
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Suppose that the received spreading code phase uncertainty is equivalent to no 
longer than L chips le time uncertainty region between u < ened and local code is 
— +^p- where T c is code chip interval Thus to avoid am imbiguity an acquisi 

tion code which is minimum of L chips long is required If L is not a pou ei of 2 choose 
a code length which is smallest powei of 2 greater than 2 Let n = lop 2 L and constant t 
a binaiv n tuples 

b, = (o) c] a], a]) j = 1 2 3 / (2 .) 

where crj" =— 1 such that 

(2 0 

These n tuples aie binary representations of number 0 through L 1 with MSB pi iced 
on the right Given these vectois the rapid acquisition sequuue is 



These sequences aie simply the tianspose of columns of a (L \ n) matrix whose rows 
aie bj In addition to high cross correlation between x and s, (he cross correlation foi 
all 6 t , i = 1 2, 3, , n aie equal and is given as 

for n odd 

(2 10 ) 

for n even 

In spite of longer integration time for each binary decision significant acquisition time 
improvement can be achieved with these codes The puce paid is reduced jamming 
resistance and increased detectibihty for the spread spectrum system Both of these 
factors are the result of additional structure (i e , decreased i mdomness) of the code 
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2 4 Coarse synchronization in direct sequence SS 
systems 

A spiercl spc c tium leccivei must generate n pseudo noise sc cjii( nee tli it is s\nchiom/( d 
with the le e ened sequence that is the corresponding chips must piee isep coincide A.s 
shown in Fig 2 3(a) the auto conclation function (ACF) of i PN sequence is i n imra 
tnangulu pulse of width two chips my misalignment cruses tin signal amplitude it 
dcmodulrtoi output to fall m accoidance with the auto conelition function A.s such 
the accurate svnchionizrtion is one of the most important functional requirement fen 
any sprerd spectrum system In an hostile environment rn requisition system with 
the following feat me is desuecl [9] 

• Bee msc successful ] iniming diumg rceiuisition coniple te h chsrbles r e omnium 
crtion system the requisition s\ stem must hwc a stiong mteifeience ie]( e tion 
c ip ability 

• The PN sequence used foi acquisition system must be piogiammable and sufh 
ciently long foi security 

• The acquisition should be lapicl so tint a jamme r must opei rte continuoush 
duiing requisition continuous opeiation reduces the uiiount of j rnimmg pow e l 
that c rn be piocluced 

Rcfei to Fig 2 1 foi a bettei understanding of coarse acquisition process m DS SS 
system If the PN sequence of receiver is synchronized with th it m ti rnsmitted sign il 


Received 



Tigure 2 4 A simple direct sequence acquisition block diagram 

then the received signrl i(t), collapses back to narrow band tiansmitted chta stre rm 
u(t) Filtering removes any remaining high frequency tains Once the signrl is 
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despiead and passed thiough the filtei u ( t ) is analyzed bv some detection method 
Howevei if the received signal and locally generated PN sequc nee aie not of the s mu 
phase the resultant output fiom multiplier and then filtei ippears as band limited 
Gaussian noise Thus to extract any intelligible data from a spread spectrum usn the 
tiansmittei and receivei must have aligned PN codes All methods of peifomung co use 
synchronization involve evaluating the correlation at \anous delay \alues M itched 
filter methods slide the signal past the PN code to pioducc the con elation R(r) foi 
all delays r mcl wait for a peak to indicate the conect code phase Other methods 
starts with a guess tq as to the conect delay evaluate R(tq) loi a possible windowed 
veision of the code see if it exceeds the thieshold(s) and upon failuie guess anotln 1 
delay estimate r 0 accordingly 

For coarse acquisition of a length L periodically repeating spreading code to i 
resolution of 1 /q th of a chip, theie aie qL possible offsets lx tween the transmit and 
receive spreading codes The goal of any coarse acquisition procedure is to leducc 
the uncertainty regarding the actual offset value as quicklj is possible subject to tin 
constraints on the available processing powei An important ic search area within co use 
acquisition is to determine the seaich strategies for searching llu remaining unceit unh 
legion [20] That is when a SS signal is received there is an issociated unceit ainti (oi 
unceitamty legion) for the locally generated PN code A tnal code phase is teste cl 
and if it is accepted code tracking (fine PN code synchroni/ ition) is invoked If tlu 
code phase undei test fails another phase is tested Search strategies are the process 
by which subsequent code phases are selected upon failure of the preceding code ph ise 
A numbei of seaich strategies aie possible depending upon the amount of a\ ail able a 
prion knowledge of received code phase 

2 4 1 Serial search fixed integration time single dwell acquisi- 
tion scheme 

Semi search acquisition consists of a seaich usually in discntc steps, of the possible 
time alignments of a receiver generated PN sequence relative to a leceived PN sequence 
[2 20] The time unceitamty region is usually quantized into i finite numbei of sc nch 
positions called cells If the search step size is denoted by s inel number of code chips 
m unceitamty region by N u then maximum numbei of cells q to be searched is gn c n 
by q = N u /s Typically s is chosen to be - so that q = 2N The cells are sen ilh 
tested until it is determined that a particular cell coi responds to the alignment of 
the two sequences to within a fraction of a chip i e , only one out of these ci e e 11s 
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iepir seats the comet code phase foi acquisition The simple model of i single rlvr 11 
semi PN ic qiusition sastem is lllustiated m Fig 2 5 which cmplcns non coheieat t\ pr 
(squue law) detectoi The fixed integration time deter toi is so aimed bcc iusr thr 



Figiuc 2 5 Single dwell serial search PN acquisition system [6] 

intcgiation pcuod (couelation peuod) of incoiiect and count cells foi this dctutoi 
leniuns const mt The icc< ivc d signal plus noise is conelited with i loc il ir phe i of flu 
PN code and then passed through a band piss (pre detection )hltei The filter output 
is then squue law envelope detected with the detectoi output mtegi ated foi fixed 
time duiation ( the dwell time ) in an mtegxate ind dump cncuit (post detrition 
mtegi ation) uid then cornpaied with a thieshold If thr thirshold is not exc titled 
the cc 11 unchi test is rejected and phase of local sequence is ldvanci d oi ictuchd to 
fest a new cdl When thr thieshold is exceeded seaich is stopped uid unfit itiou 
is clone Once the puticulai cell passrs the venfication tc si control is handed rnn 
to ti aclung circiuti} foi fine: alignment of the two codes [G 9] One of the obi nms 
disadvantages of this rouse acquisition technique is that when 1 ugc unceit aintv legions 
uc involved a lot of tune may be spent unnecessanh intr gi ding the mcon.ec t corlr 
phase cells Flowevci it is most secui e acquisition techmcjut m hostile environment 
due to its extiemely stiong mtciference rejection capability [9 20] Moieovei it is lit st 
suited foi low SNR spiead spectium svstems when compaurl with othei acquisition 
techniques [15] 
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2 4 2 Multiple dwell acquisition scheme 

The multiple dwell detection scheme was designed to leduce the time spent in intcgi it 
mg the meoiiect code phase cells [22] A diagram of multiple dwell detector is shown 
m Fig 2 6 as taken fiom [6] 



Figure 2 6 Multiple (N) dwell PN acquisition system [6] 

The s'n stem begins by selecting a local code phase for e-s dilation multiply it w ith 
received signal plus noise The envelope of lesultant product signal is formed m i 
square law envelope detector The envelope is then applied to c uh of the ATmtegritois 
simultaneously The integration (“dwell”) time of detectors uc oideied such that 

Tdl < T d 2 < r d3 < < T dN (2 11) 

Thus all of the N integrate and dump circuits starts integration at the same time, c xch 
one dumping, however at a later time instant in view of (2 11) The seaich updite 
algorithm woiks as follows 

The output of i th integrate and dump is sampled and comp ued to a threshold only 
if all of the previous % — 1 integrate and dump outputs have previously exceeded the n 
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respective thresholds Otherwise the first integrate and dump output to fall below its 
threshold causes the local code phase update and search the iu \t cell thereby resetting 
all of the integrate and dump circuits The maximum turn to seaich a giaen coll is 
TdN while minimum time is t^i Thus most of the cells can !h dismissed after a dwell 
tune l < IV which implies significant savings in the acquisition time as opposed 
to simple single dwell sennl search algorithm where each all is examined foi i tnnr 
equivalent to t,^ 

2 4 3 Sequential estimation method 

This scheme was first introduced by Ward [25, 26] and is bas< d on a sequential estirna 
tion of shift legistei states of PN generatoi In particular tin RA.SE ( Rapid \cqui 
sition by Sequential Estimation) of Fig 2 7 makes its best estimate of fust n leant d 
PN code chips ( n is the number of shift register stages in the loc il code geneiatoi) uid 
loads the rccener sequence generator with that estimate thus defining a paiticulu 
initial condition (starting state) from which the generator begins its operation Since 
a PN sequence has the property that next combination of register states depends onh 
on the present combination of states If all n input chips ait c stmiated w ith suflkit nt 
accuracy all the following states can be predicted based on the knowledge of onh tins 
initial condition Foi this the local PN code generator is periodically loaded with new 
estimated n chips at a rate determined by examination penocl generatoi Win n to 
stop the loading is determined by a threshold crossing of m in lock detector The 
test statistics of the in lock detecter is the cross correlation of input code w ith that 
produced by local PN generator whose shift register contents < ouespond to the piesent 
estimate of the input state Once it is determined that the n d< tccted chips pioy ide the 
coiicct staiting positions lock detector threshold crossing inhibits fuithei loading of 
PN code geneiatoi This register is then closed within tlu PN tiaclung loop which 
is responsible foi maintaining code phase from that time onw uds This scheme out 
performs the single dyvell serial search method for modeiate ( lOdB) to high (+10dB) 
input SNR [26] However its performance is inferior to send sc uch at yen low SNR 
(i e below 10 dB) [6] Despite its Rapid acquisition capability this scheme is highly 
vulnerable to noise and interfering signals [30] 

2 4 4 Sequential probability ratio testing 

Sequential Probability Ratio Testing (SPRT) is a direct seqmua acquisition method 
which is in some ways quite similar to the multiple dwell detc ( ( oi Initially tm estig itc d 


Chapter 2 


21 



Figure 2 7 A rapid acquisition by sequential estimation (RASE) technique 

by Wald [2 6] it has received considerable attention in recent u \rs [27] Fig 2 8 depicts 
a block diagram of an SPRT acquisition system The difference between SPRT uicl 
other acquisition scheme lies m the decision circuitry SPRT uses a two level (A uid 
B) threshold device whereby three possible events occui The se being 

• If the paitial con elation is less than threshold A, declue code phase as me one c 
t tiv a new phase and leset the integrator 

• If the partial correlation is gi eater than Threshold B(> 4) then dechic cone cl 
code present 

• If the partial correlation falls between Threshold A and B then retain the conect 
code phase and integrate a little longer 

The dual threshold system implements a screening test Test lc suits which cleail} show 
success or failure of acquisition allow a quick decision Those which are inconclusive 
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Figuie 2 8 Sequential probability testing block diagram 

aie used to initiate further investigation m this case longei mtegiation 

As can be seen from Fig 2 8 the received sequence is conelated with the localh 
generated \ersion of the PN code, product signal is square lav envelope detected and 
a bias le\el is subtracted to remove any biasing effects fioui the squanng of noise 
The residual signal is then integrated and sampled at inter \ ils (oi multiples) of dnp 
period T c The sampled signal is then tested against afore m< ntioned thresholds uid 
appropriate decision is subsequently made SPRT operates on the principal design 
philosophy of discarding false synchiomzation propositions quickly and accepting tine 
ones quickly 

The si stem descnbed above has many negative aspects Fnstly there aie moie tli in 
one threshold and a bias level from optimization point of view [6] Secondly the ts pu il 
SPRT si stem may suffer fiom the significant problem of sign U remaining between 
threshold A and B indefinitely The Truncated Sequential Probability Ratio Testing 
(TSPRT) was introduced by [27] to address this shortcoming TSPRT essentially sc ts 
a limit on the longest correlation period and if this limit is leached a hard decision 
is made whether the correct code phase is acquired or not A significant chawbacl of 
(T)SPRT is that it is difficult to analyze due to evaluation of paitial correlations uid 
the variations in integration period Consequently mean uid vuiance of acquisition 
time is best obtained via simulation [6] 
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2 4 5 Matched filter approach 

All of the pie\ious discussions of DS acquisition techniques mcorpoiate a measun of 
coirelation which was produced by an active coirelation of the received signal with i 
locally generated PN lefeience code The term active condition is used w ith multiply 
and integrate type of correlation/detector structure wheie n (( i\ed signal is multiplu d 
with the local PN code energy detected tested against a thnshold(s) and declaicd as 
a conect oi mcoirect code phase Moieovei the local PN gtneiator is continuoush 
running theiefoie foi each successrve threshold test a completely new set of 1/ = rjf- 
chips of received signal is used foi each successive threshold test On the othei hind 
matched filter based acquisition system operates on the pnnoph of passu e conelation 
where the input (received signal) continuously slides past i stationary replici of PN 
sequence when there is a match (i c the two are m near synchronism) the threshold 
is exceeded local PN generatoi is enabled and despreading cm commence If the c ode s 
do not match the threshold is not exceeded and PN generrtion is not started The 
mam advantage of matched filtei approach over the previous methods is summan/cd 
m [4 6] 


The reduced synchronization time is explained by noting that the MF 
continuously correlates over a time interval That is tin MF merely waits 
for the preselected waveform to be received 

The base band matched filter acquisition system is shown in Fig 2 9(a) A tians’ursal 
matched filtei peifoims many more possible code phase conclations in an internal 
than any of the preceding methods due to its inherent panllel stiuctuie which does 
not requne a new' observation interval of received signal foi c \ch integi ttiou peuod 
As a lesult the search rate of matched filter acquisition method is considei rbh higher 
than piecedmg methods [1 6 21] The matched filters aie fibnc rted m the following 
forms 


• Charge Couple Device (CCD) 

• Surface Acoustic Wave (SAW) Convolvers 

• Discrete Time Conelators 

The last is becoming increasingly popular due to the availabilitv of high speed Digit il 
Signal Processors (DSPs) and flexible Field Programmable C ate Arrays (FPGAs) A 
digital correlatoi implementation of matched filter is shown m Fig 2 9(b) Using this 
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Figure 2 9 (a) Low pass version of matched filter acquisition (b) A digital 
correlator implementation of matched filter 

implementation of matched filtei m Fig 2 9(a) results m the discrete time leision of 
the matched filter acquisition system The correlation output is then squaied, summed 
up and finally thieshold tested to decide whether the code is u qunod oi not Matched 
filter technique is prefeired when short PN codes are used 


2 5 Fine synchronization m direct sequence SS sys- 
tems 

Once the coaise bjnchionization of PN codes is achieved code Tiaekmg is activitid 
The purpose of code tnclung (fine PN code alignment) is to nn\imize the aligimunt 
of locally generated PN code to that of received incoming PN code and in doing 
so to maximize the post despreading SNR In general cock flicking loops maintun 
synchronization of leceiver s replica of the spreading code by using the two correlitois 
called an eorly correlator and a late correlator An early coiulitoi uses a code iefeienu 
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waveform that in time by some fraction of a chip with iespcc t to acquired Pis code 
A late correlator uses a code reference waveform that is delucd by some fraction of 
a chip The difference between the two correlations is then used to form error signal 
proportional to the code misalignment and uses this enoi sign il to shift the two cocks 
mto closei proximity There are two basic architectures which implement the required 
fine closed loop functions [1 4 6 30] the Delay Lock Loop (DLL) and the Tau Ditlioi 
Loop (TDL) 

2 5 1 The delay locked loop 

A delay lock loop corrects any timing offset of the local code generator by feeding 
back enoi \oltxges to the local code generator Considei the simplified block dngi un 
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[3] as shown in Fig 2 10(a) Timing eiror voltages aie obt uncd by using th< Io< il 
code to generate two offset sequences (advanced and dplai < d) These offset < ode s 
aie separately multiplied with the input m the two channels of delay locked loop 
The multiplied outputs ue then band pass filtered and emclope detected md non 
coherently combined to produce a difference signal e(r) wlnu e(r) = I 2 (t) — 1 1 (t) 
This difference (or correction) signal is filtered (dampened) 111 the loop filter and is 
used to derne a voltage controlled oscillator (VCO) The \ CO is a clock which d 1 1 \ f s 
the local PN geneiatoi fastei when its clock is lagging in phase m comparison with flu 
incoming sequence and vice versa The rigoious treatment of DLL is given in [1 4] md 
a simplified veision in [3] The normalized error signal as obt lined in [3] is given In 

E(r A) a R\ (t - - R* (r - (2 12) 

Where R x and i ?2 are the auto correlation functions in the < rrlv and lxte bimch of 
DLL Parameter r is delay estimator error , defined as r = T/ — T d T<i is the clelax 
estimate and A = separation between early and late codes Lndei the assumption of 
identical blanches R x = R 2 = R 

E(t A) = & (r ~ f) ~ R 2 (t ~ f ) (2 13) 

The auto correlation function as given bj (2 12) and (2 13) ch pends on the basic chip 
waveform (j) c (t) For a rectangular waveform of unity amplitude in the interval [0 T,] 
the error signal is shown in Fig 2 10(b) By observation to the Fig 2 10(b) it can Ik 
concluded that regai dless of whethei A is positive oi negatnc a compensating <noi 
voltage is generated which will correct the chip timing difleience However it should 
also be noted that DLL is difficult to construct since the advance and delay aim of 
stiuctuie must be peifectly balanced requiring identical nn\c is bandpass filteis uul 
envelope detectors 

2 5 2 The Tau- Dither loop 

To overcome the balancing problem of delay locked loop the Tau Dither Loop (TDL) 
is suggested TDL operates on the principle of time sharing a single correlator bi c uh 
and late codes le the incoming signal is correlated alter nntehj with the eaih md 
late codes [3, 2, 30] The Tau Dither Loop is shown Fig 2 11 [3] This airangenu at 
generates the same code offsets as did the previous delaj locked loop, howeiei onh 
one multiplication branch is used at a time periodically alternating between adi anted 
and delayed version of the code eveiy T p seconds wheie T) is a fraction of \ chip 
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Figuie 2 11 A Tau Dither Loop for PN code tracking 


period (usually taken as i) The alternation is achieved by means of a gating signal 
q(t ) The gating signal also used to provide the resultant erioi signal before smoothing 
(averaging) in the loop filter 

The nonnalized error functron error function as given bt [3] is 


E{t t A) = R 2 



This enoi function consist of two parts the first part is famili n enoi function of DLL 
as plotted m Fig 2 10(b) and the second part rs sum of the squires of the early mcl 
late auto correlation functions multiplied by the dithei swift lung function q(t) which 
causes a net reduction in the error correction voltage 
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The Single Dwell Serial Search PN 
Acquisition System in Low SNR 
CDMA Systems 

3 1 Introduction 

In a tjpical spread spectrum system the requirement of laige processing gam is ic 
comphshed by using high chip rate PN codes Coupled with this is the requnement 
of a high detector input bandwidth which in turn results in low input SNR Channel 
distortions such as multi path fading causes the received signal to be random The 
multiple access application of spread spectrum system such as land based mobile CD 
MA. netvoil (e g IS 95) or Low Earth Orbit (LEO) satellite s's steins (e g Globestar 
etc ) suffer fiom another phenomenon called Doppler effect which is due to rehtnc 
motion between transmitter and/or receivers It has been shown m [1, 15] that SNR 
at the conelatoi output will be drastically low if theie is considerable Doppler Low 
link power budget is another reason for the recened SNR to be low In general m 
any CDMA network the excessive number of users approaching the system capnciti 
limit causes the received signal to be characterized by low SNR Since this SNR ilso 
happens to be a key parameter influencing the decisions regarding the synchronization 
therefore code acquisition process will also suffer heavily The t\ pical operational sc e 
nano m terms of input SNR at receiver is well below -10 dB [6] however, above facl ois 
further i educes it considerably causing severe degradation m system performance 

In this chapter the serial search code acquisition scheme of [15] has been moclifn cl 
foi extending the low SNR operation range 
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The code acquisition scheme of [15] consists of a multiphei band pass filter recntu 
lation loop accumulatoi a decision circuitry acquisition and t lacking contiol cirnntir 
and a local PN code generator The scheme has shown good ipsults in the SNR langc 
— 10 dBto + 10 dB fiom acquisition parameters point of \iew i c an impio\ement of 
3 times m mean acquisition times and 5 times in the vaiiancc of acquisition time h is 
been obtained compaied to standard serial search approach of Fig 2 5 Compai iblc 
performance impio\ement is also obtainable by sequential estimation methods such is 
RASE [2d] 

Howevei below —10 dB SNR the system performance needs to be le investigated 
since improvement is nominal (1 5 times) and for this purpose two major modifications 
m the scheme of [15] have been introduced with due consideration to hardware oui 
head The first one being the addition of post detection integration (PDI) blocl 1 xncl 
the second one being adaptive threshold operation for achieving the constant false alar 
m probability using cell averaging constant false alarm rate (CA CFAR) approach 2 
m contrast to the ADTLC algorithm [17] used for the same pm pose in [15] 

The attractive feature of the scheme presented in [15] nunely the recirculation 
loop is retained in this woik and after the above said modific it ions the modified scual 
seaich code acquisition scheme is piesented m section 3 5 The analysis of modified 
scheme is then earned out m subsequent sections from detection and acquisition point 
of view 


3 2 Recirculation principle 

The recirculation loop operates on the principle of delayed addition snnilai to tint ol 
a Fust In First Out (FIFO) accumulator queue The samples of decision vanable 
are fed to decision piocessoi and the accumulator simultaneously The delai piouded 
by the accumulator can be chosen according to the desned ipphcation As shown in 
the Fig 3 1(a) the delay m accumulator is kept T p = qr, d win u T v is the period of the 
con elation pulses at the output of sliding correlator as shown m Fig 3 1(b) q is the 
number of cells in the uncertainty region of code phase t/ is the integration (dwell) 
time Therefore, lecnculation loop accumulates the signal samples in each cell obtnued 
during every new letrace of code delay uncertainty region Th< addition of recirculation 
loop in standaid single dwell serial search approach of Fig 2 5 will impiore the SNR 

1 The necessity of PDI is discussed in section 3 5 

2 CA CFAR method is discussed in section 3 3 and 3 4 
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(b) 

Figure 3 1 (a) Accumulator recirculation loop (b) correlation pulses at the 
output of correlator 


of decision variable which will lead to better probability of s\ nchronous cell detection 
and improved acquisition time performance Simulation of ucnculation in base band 
foi L = 127, q = 2L is shown in Fig 3 2 


3 3 Adaptive CFAR approach 

The basic cell avuagmg constant false alarm rate (CA CFAR) piocessor is lllusti it t d 
m Tig 3 3 The essence of CFAR is to compare the decision st itistics with an adiptiae 
threshold, obtained from the estimate of noise oi interference power aftei obseiung 
fixed nurnbei of noise or interference samples This estim itc may be obtained In 
accumulating the noise or interfering signals and adjusting the threshold in record nice 
with the measured noise oi interfering signal so as to provide constant false alum 
probability As shown in Fig 3 3 the samples of decision \ uiables are fed to a shift 
register and it is assumed that only one memory element of shift legister contains signal 
plus noise sample while remaining elements contain only noise s amples An estimate (X 
) from these noise samples is then formed and multiplied by a thieshold coefficient (T) 
T keeps on changing according to estimated noise (X) so as to foim a threshold level 
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Figure 3 2 Recirculation loop output with q = 2L , L — 127 and recirculation 
delay T p = 254 LT 

rj — T\ agunbt w hich the samples of decision variable are complied and a decision is 
made whether signal is present or not More details about CF4.R appioach are „i'en 
in section 3 6 2 

3 4 Adaptive CFAR vs fixed threshold approach 

The thieshold level m the decision device of anv commumc ition receivei is chosen so 
as to achieve a desued false alarm probability The false alum rate is quite sensitive 
to the thieshold level and minor variations such as changes m the noise at the input of 
receiver multiple access interference shadowing and fading m the leceived signal ninv 
affect the receiver s detection performance drasticallv The tv pu il mobile environnu nt 
is further characterized bv the relativ e motion between the ti uisnutter ancl/oi lecuvc i 
which adds a new dimension to alreadv comple\ code acquisition process In the 
fixed threshold detection methods the threshold level is fixed for a given false alum 
probability no matter how the associated environment of it reived signal is v uvm 0 
and is changed onlv when a new value of false alarm probibilitv is desired It will 
quite often lead to too mam false alarm which wall reflect in an increased acquisition 
time resulting the poor acquisition time performance [17 18] On the other hand 
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Figuie 3 3 A CA CFAR based decision processor 

in adaptive thieshold detection methods the threshold level is continuously upditcd 
such that it piovides desned false alarm probability by foiinmg an estimate of the 
noise or interference m the received signal [18] Therefore the problem of automatic 
threshold level control is fundamental for efficient operation of a DS SS sj stem and 
forms a key issue m decision making about the correct code phase (acquisition) One 
such approach for adapting the threshold level is Automatic decision threshold lc \ el 
control (ADTLC) [17] which comes at the cost of fairly merexsed hardwaie complexity 
coupled with the fact that it is SNR dependent and therefou not well suited foi low 
SNR environment Less demanding hardware is needed in the cell average const ml 
false alarm rate (CA CFAR) methods [18 31] and they aie also independent of iecei\ c d 
signal s SNR 


3 5 Improved Serial search hypothesis testing mod- 
el for low SNR CDMA systems 

The recirculation loop of Fig 3 1a post detection integration block and the ad ipt ivc 
decision piocessor of Fig 3 3 has been added in the standard st rial search appioich of 
Fig 2 5 

The resulting modified serial search code acquisition system is shown in Fig 3 4 
and is expected to have better performance in terms of detection probabilities and code 
acquisition panmeteis as compared to the standard senal so ik li approach and tin one 









Chaptei 3 


33 


piesented in [15] due to the following leasons 
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Figure 3 4 Code acquisition loop employing recirculation principle and adap 
tive decision processor 


• Post detection integration (PDI) block is thought to lx necessan foi low SNR 

operations because the key paiametei influencing the conect decision is the S 
NR at the multiplier output It has been shown to lx piopoitional to squuc 
of code conelation function multiplied by a factor of due to cunoi 

Dopplei wheie u > d is cainer fiequency shift in ladian due to Dopplei effect md 
Td is integration or dwell time [1 15] The SNR at multiplier output will be 
drastically low if theie is considerable Doppler Even under the assumption of 
no Doppler testing the PN sequences with half chip timing erioi would c uisc 
squaie of code correlation function to be quite low and will result m seiei il clB 
loss m the performance of acquisition scheme a hardh affordable factoi in low 
SNR environments PDI forms unbiased estimate ovei the multiple sample s of 
same hypothesis which in turn impioves the SNR of decision \ an able then b\ 
justifying the necessity of post detection integration 

• The adaptn e CFAR decision processor being mdepende lit of SNR, has been use d 
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to pio\ide a constant false alarm probability (Pn) dining eich reti ice rule of 
seal ch procedure 

The unpio\f ment m SNR dm to recirculation loop and post d< t< ction mtegi rtion done, 
with constant false alaim probability will increase the overall detection probability m 
each of the uncertainly legion cell during evei} letrace cvclr 

The serial search process for code acquisition comprise s of multiplication of recent d 
direct sequence signal with the local PN code sequence meisuiement of eneigv m the 
pioduct signal and then comparison with a threshold in the decision processoi The 
samples of decision variable Z x are also fed to the accumulitoi lcciicul ition loop If 
the two sequences are not aligned band pass hltei blocks most of the power in pioelue t 
signal and tlueshold is not exceeded which results m the lejeetion of code under te st 
The phase of local sequences is then retarded or advanced foi te sting a new code ph tsc 
Thus c iclr code phase correspond to the two hypotheses when the two sequent ( s ue 
lie ally aligned the pioduct becomes a narrow band signal tli it passes thiough the 
filter with little loss of power anel threshold is exceeded indie itmg that code has lx e n 
acquired anel search is stopped coirect code verification is clone and the ti iiling 
system is activated This procedure is repeated for every < ode phase (cell) m flu 
uncertainty legion If cluimg the first pass the threshold is not exceeded then duiui^ 
the next retiace cycle of uncertainty region samples of decision variable are added 
to the corresponding samples of previous search cvcle thereby improving the SNR of 
decision v enable which will in pimciple impiove the acquisition time performance 

3 5 1 Probability density function of decision variable 

The complete statistical description of decision variable Z of Fig3 4 is given by its 
probability density function First the simplest case is considciecl m which no post 
detc ction integration is employed Next post detection mtc gi ation case is employe el 
In both the above cases the received signal is assumed to lx unfaded Finally tlx 
behaviour of decision variable under Nakagami m faded lece iv e d signal is conside k d 

3 5 11 Unfaded signal m AWGN No post detection integration case 

Let the input data waveform is denoted by d(t) The transmittc cl signal aftei spreading 
by PN code wavefoim c(t) can be expressed as 


(3 1) 


s(t ) = c.(t)d(t)VPPcos(ujQt + a/') 
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Wheie ixio is cmiei frequency in iadians and P is the pova m tiansmitted cinui 
signal The locened signal in Fig 3 4 is r(t) = s(t) + n(t) vlic k n(t) is assumed to b< 
idditne white gaussian noise The local PN code geneiatoi ol Fig 3 1 is assumed to 
be at an offset of ±r fiom the leceived code so that local PN code is lepiesentcd is 
c(t ± r) In the woist case scenuio r may be as large as the pt nod erf PN (odt 1 he 
lecened signal altei conelation and band pass filtering can lx expressed as 


t(1) = [r{t)c(t±T)} BPF (j 2) 

wheic [ ]bpf represents band pass filtering operation The input to squne law de tcctoi 
then can be cxpiessed as 

r ( t ) — V2 ~P p(f)d(f)cos(o/ c ,f + 1 /’) + V2n c (t)cos(tx> 0 t + ip) — \f~2n + ij ) ( > 3) 


where pit) is conelation function of PN sequence and is define d is [17] 

^ _ | 1 if the two PN codes are synchronized 1 e signil lit s on the conelation emu 
0 othciwisc 


Letting 1 = p{t)d{t)\fP and lccogmzmg A as ims signal amplitude after coirehtum 
and filtering is performed P is the power in the leceivccl sign ll iftci conelation aid 
band pass filtering 

Now defining the hypothesis Hi as the case of signal pnsait it (the cell bung 
searched corresponds to a sample value on PN correlation cuive of Fig 1 1(b) and 
p(t) — 1) and hypothesis H 0 as the case of signal absent 1 < (the cell being searched 
docs not correspond to a sample value on PN correlation curve of Fig 3 1(b) and p{t) = 
0 ; 

The input to sqtnic liw envelope detectoi of serial seaich hypothesis deuce of Fig 3 4 
can be expressed 111 the foim 


%(t) = V2 Acos(uj 0 t + ip) + V2n c (t)e.os(u) 0 t + ip) — \/2 ?i {t)sm(oj„t + iji) 


(3 5) 


— s/2R(t)cos(uJ 0 t A - ip Ac @(t)) ( 3 6 ) 

Where 


R(t) = vW”c(<)) 2 +»*(#) «(() = fan' (3 7) 

n c (t) and v s (t) aic band limited independent low-pass mo mean Gaussnn noise 
processes with vanance a 2 = ^ N 0 is single sided noise spectral density and B is 
the noise bandwidth of predetcction band pass filter 


The output of squ up lav pm elope detectoi m Fig 3 4 m i< spouse to the input / (t) 
of (3 6) is (ignomig the second harmonics of camei) 

y(t) = ^(t) = R 2 (t) = {A + n (t)Y + n (t) ( 8) 

Letting 


4 + v (t) = w(t ) ~ Af{ 4 <j‘ > ) 
n $ (t) = v(t ) ~ J\f(0 a 2 ) 


( 0) 


the distnbution of j(t) under hypothesis H% is found out in ippnidix B 1 


i e 


Pi(y) 


= I 2T ex B (- 2^ + 7,) h (7?7) *01 <J > 0 


oth( nvisc 


(3 10) 


Which is non central chi squared pd/chaiactenzmg the squ ut 1 m output m all s< nth 
cells I hat contain the signal to be acquired embedded in YY\ ( N is signal to noise 
i atio m the i th retiace cvcle of recirculation loop defined as 7 = /70 70 is piecletation 
signal to noise lutio given by 


a A 2 _ A 2 
~mT 0 ~ 2 a* 


(3 11) 


Iq{ ) is modified Bessel function of first land and zeroth 01 da 


The pdf of v(t) unclei hypothesis H 0 can be obtained hs substituting 4 = 0 in (3 1 0) 
and iecogni7ing the fact that I 0 (0) = 1 


po{y) = 



fol (/ > 0 
otheiwisc 


(3 12) 


This is antral chi squared pd/ characterizing the square lav output m all seaicli c< 11s 
th at cont am C aussian noise only Equation (3 10) and (3 12) u< same as those giic u in 
[6] except the fact that definition of 7 is diffeient heie to mount foi the icarculatiou 
loop 


3 5 12 Unfaded signal m AWGN Post detection integration case 

If band pass filtei of Tig 3 4 has sharp cutoff skirts and v(l) is sampled at mten als 
T — then the samples (y x y 2 dN B ) will be approximately independent [i 6] 
Therefoie integrate and dump (post detection integration) output can be approxuu at 
eel by a summation ovei these sampled values, namely 

z = - y{t)dt » — 2 ^ % 

Td 4o N b k=Q 


(3 13) 
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N b = j ; = Bt 1 (ill) 

wheie Td is dwell time Th pdf of normalized decision \amble Z = (denied m 
appendix B 2) undei the two hipotheses H 0 md H x lespectneh can Ik be e\pie sse cl 

' ) fm z >o 

0 olheiwis< 

and 


p (Z ) 


(3 1 j) 


Pi (Z ) 




exp(-Z - Ngril^i&y/NTrZ ) foi Z >0 

otheiwise 


(3 10) 


3 5 13 Nakagami m faded signal m AWGN 


Often m the tcnestml propagation signal Aimes at the Hienei tlnough drift k ill 
p itlis with lelitive clelais and l indom phases Statistic alh the usult nit sign il tin 
phtuclc is modeled bv the Rayleigh oi Rician di&tnbution while pln&e is assumed to Ik 
umfoimly distnbutcd m the mtervxl [0 2ix] Howevei as poinle d out m [32] moie com 
plete md accurate statistical description of urban lacho multipith fading channels is 
provided by Nakignm m distribution as compared to Rayleigh oi Rician distribution 
Also the Nakignm m is a two parameter chstiibution tint pi exudes more flexibilite 
and accuracy m matching the observed signal statistics as compared to the othu two 
chstiibution Foi fading parameter m < 1 the Nakagami m distribution is used to 
model fichng channel conditions that aie moie se\ ere than the Rayleigh distribution 
( m = 1 ) and for m > 1 it is used to model the fading climnel conditions less se 
ve re than Rxxleigh chstiibution It has been xs&umed m the following discussion th it 
Nxkxgxim m fading is slow enough that amplitude as well is phase lemam const uit 
ovei the one mtcgialion ( dwelt ) time (r^) but fast enough th it successive integration 
dm at ion in essentixlh independent The Nakagami m fading chstiibution is gru n In 


[33] 


p(U = 


2 


m 


r(m) Vo 

wheie m is fading paiametei and is defined as 

O 2 


i 2m l exp 


-m j 

IT 


m = 


(3 17) 


(3 IS) 


E [a 2 - O] 2 

and O = E(x 2 ) and T(m) is gamma function In Nakagami m fading enviionme nt 
the received SNR will be a random variable with chstiibution obtained by a change of 
vmable m (3 17) as 


P{h) 


% 


,771 — 1 


m\ 

r(m) VO J 


m 


exp 


myA 
' 0 J 


(3 19) 
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Even m scveit fading environments p 0 (Z*) as given by (3 15) will lemam unchanged as 
it does not depends upon the fading signal power Howec ei p\(Z ) as gi\en b\ (3 16) 
will now become conditional pdf 1 e P\{Z /* y ) Theiefore pdf of decision variable Z 
undei Nakagami m fading will then be obtained by avenging the pelf (3 16) ovei (3 19) 
l e 

r 00 

Pi{Z)= Pi(Z /y t )p(y t )d ll (3 20) 

J 0 

This integi vtion has been earned out in the appendix B 3 md the final lesult c m In 
expressed as 


{Z ) Nb ~ 1 exp[-Z*) ( 

_ r(m)(l + u) m (N B - 1)1 lFl ( m Ne 


uZ 
1 + u 


In pothesis Hi (3 21) 


wheie ir’i(a b x) is the confluent hypergeometnc function given as iFi(a b i) = 
S£To dn ^ ^ ^ ~ 2 [14] u is defined as » = anc | Q = E[ , ) 

Since u is a function of average SNR in the i th retrace ejele so is the Pi[Z ) 


3 5 14 Rayleigh faded signal in AWGN 


The Rayleigh fading distribution is a special case of Nakagami m distribution [33] it 
can be obtained by putting m = 1 m the latter distribution Therefore for pdf of 
normali7ed decision vanable Z* under Rayleigh fading is given by 


*(* > - h n b 


1 + u 


1 + v 


In pothesis Hi 


or altcrnativelv by 


Pl(z) = (£&(z£ if 


(ft) 


1 + u 


r=0 (Nb + 1 ~ l) 1 


(3 22) 


(3 23) 


The lmpoitant case of no post detection integration case undei hypothesis Hi can be 
obtained by putting Ng = 1 in (3 23) The final pdf is given is 


Pi(Z‘) = 


1 4" 


exp 


Z ' 
1 + 


(3 24) 


With Z* now given as Z = This result is similar to one obtained m [1] except 
the change of notation 


3 6 Evaluation of detection and false alarm proba- 
bility m terms of PN acquisition parameters 

The detection probability, P D is a salient factor in acquisition time and is defined 
as probability that the code alignment is declared when the code is aligned If Pp is 
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low it is possible that multiple sweeps of all possible code phases max be necessaiv 
before acquisition occuis Probability of false alarm Pfa /*> the probability that the 
correlation product of the received and reference code is sufficn rit to declare acquisition 
when m fact the correct code phase has not been acquired The event of false ilum 
results in futile attempt to despiead over another integration period Both probability 
of false alarm and detection aie parameters of code acquisition and acquisition tunc 
Tacq The Pq and Pi t aie also functions of the mtegiation pc nocl ( dwell firm ) 
SNR and threshold voltage The ideal scenano is to ha\c high Pp and love Pi i 
which is not always achievable The variables which affect the process of determining 
these probabilities aie intei lelated so that a simple change m the mtegiation tunc ox 
threshold voltage affects both of these probabilities A high Pp and low Pf a maph i 
long mtegiation period and consequently long acquisition time 

3 6 1 Unfaded signal m AWGN Fixed Threshold method 

The probability of false alarm Pfa f or unfaded signal tn AWGN and post detection 
case can be obtained by integrating the likelihood pdf P 0 (Z ) xs given In (3 15) 

/»00 

Pfa = / Po(Z )dZ* (3 25) 

Jj] 

The integration (3 25) is earned out m the appendix C 1 and result can be expressed 

aS Nb ~ 1 n p 

Pfa = exp(-r ,*) E A (3 26) 

p - 0 

where if is noimalized threshold and is defined as rf = ^ 

The probability of false alarm under no post detection migration case can be ob 
tamed b) substituting N B = 1 m equation (3 26 ) i e 

P rA = exp(-7 f) ( 3 27 ) 

The probability of detection in the fh retrace cycle of imuuhtion loop P D t m 
be obtained by 

P Dl = [ V\{Z )dZ (3 28) 

h 

The complete derivation is carried out in appendix C 1 and the final ixpiession foi 
detection piobability is given as 


(NbJiY Nb ^- in ) p 

Pdi = exp(~7 iNb -V )2s A 2^ „i 

r—0 1 P~0 1 


(3 29) 
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The detection piobability foi no post detection integration c tsc mav be obtained bi 
putting N b = 1 m equation (3 29) i e 


Pd = exp(- 


~ T ^ ( y 

■ 1 '-* )EttE“ r- 


r — 0 ^ 


(3 jO) 


p=0 


which if de&ned can be expiessed in Marcum Q function hum The lesult w ill th( n 
be same as the one used m [15] 


3 6 2 Unfaded signal in AWGN Adaptive Threshold method 


In the CA CFAR approach (Fig 3 3) for adaptive threshold s< ttmg the samples of Z 
are serially sent into a shift registei of length M +1 [31] Each of the M reference cell is 
assumed to contain Caussian noise while only one cell corresponding to Zq is assumed 
to contain the signal plus noise Non coherent integration of Nr pulses follow the 
square law envelope detection for comparison with the adaptive threshold Integrated 

= , V/ 

noise povei estimate fiom L = 1\IN B samples is then obtained b\ \ = E l= _ w Z The 
threshold TX varies according to variations in the noise and the threshold coefficient 
T is used as a scale factor to achieve the desired false alaim piob ability foi a given shift 
register of window size M Under the assumption that suiioundmg M cells cont uns 
independent Gaussian noise with the same variance, the pdf of \ is obtained fiom 
(3 15)as 

p(\) = 7pZyexp(-A) (3 31) 

The detection probability foi a given value of T is obtained b\ integrating the detection 
piobability conditioned on a threshold rf — TX ovei p(Y) [31] as 


/ do rco 

P Dl = p(X) Pl (Z*)dZ*d\ (3 32) 

Jo Jt\ 

while threshold coefficient T is determined from a designated false alaim piob abilita 
given as 

roo r oo 

PrA — / P( A)/ PQ (Z)dZ*d\ (3 33) 

Jo Jt\ 

Substituting the (3 31) (3 15) m (3 33) (3 31) (3 16) in (3 32) and carrying out the 
required integrations (see appendix C 2) the final expressions uo 


Ppa = 


1 


(1 + T) mn b 


N b - 1 

E 


p=0 


MN ft + p — 1 

V 


T 1 

(i + Typ 


(3 34) 


Pdi = 


expt-^jVe) A {NruT ( mn b +P ~ 1 ) ( T V 

(1+T) MN ‘~ 0 r 1 ~o \ V )\1 + Tj 


(3 35) 
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wheie notation 


a 

b 


is used to denote vhen a and 1) ue positive mtegei uid 


also loi icpiescntmg f(/ - +i)rff b+i.) when a and b are positm non integer 

The piobability of false alarm for no post detection integration case may be obtained 
by substituting Nb = 1 m (3 34) i e 


Pfa — 


1 

(1 + T) A/ 


(3 36 ) 


The detection probability for no post detection integration rase may be obtained In 
simply putting Nb = lm equation (3 35) i e 


exp(- 7 ,) “ 7 [ ’ / M +p - 1 \ / T \ v 

"a +Trto r 'h\ p 


(3 37 ) 


3 6 3 Nakagami-m faded signal in AWGN Adative threshold 
method 

Substituting (3 21) and (3 31) in (3 32) and canymg out tin uquired integration (sec 
appendix C3 ) the expression for detection probability m the i th retrace of scnch 
procedure is found out to be 

m + r - 1 \ / u y Nb *ZT 1 ( L + p- l\ / T \ p 
r ) Vl + n/ \ p ) Vl 

(3 38 ) 

The detect lou probability m (3 38) is a function of leceived SNR in the letxace c\ c h 
of scml search pioceduie (thiough parametei u) fading pamneter m numbei of post 
detection integration Nb and the number of reference cells M used in the adapt ru 
CTAR tluesholding method 




Di 


(1 + i/) m (l +T) J 






3 7 Signal flow graph based analysis 

As pointed out earliei, the time to acquire the code {acquisition time T acq ) of an 
acquisition process is characterized by first two moments of ac quisition time i e me an 
(T acq ) and variance {a] cq ) The two parameters are function of P FA and P Dl vlnch 
are already derived m previous sections As a next step we Mill now determine these 
two moments In the hteiature, there exist three approaches to analyze a code sc uch 
process fiom acquisition time point of view 
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• Signal flow graph based appioach [2 6 16] 

• Cncular slate diagram approach [19 20] 

• Puic probabilistic appioach [4 10] 

Foimei two approaches are z transform based and rely on st nch process descuption 
by a suitable flow graph or cnculai state diagram Next tiansfer function of gi iph 
(or state diagram )is deteimmecl Finally this transfer function is used to detcinimc 
the moments of acquisition lime These approaches are quite yersatile elegant md 
simple to use as opposed to puie piobabilistic appioach which is time domain bisirl 
heuristic in uature and lacks mathematical elegance The generalization of found 
two approaches to an arbitrary search process is stiaightfoiw arcl as they complete h 
describe the anatomy of search procedure while latter does not provide much insight 
into actual search piocess and its generalization to an arbitim search strategy is in 
involved task 

In this work modified senal search system as shown m Fig 3 4 will be mih/ed 
using the signal flow graph based approach The analysis piocedure and subsequt lit 
expressions obtained arc general enough in the sense that they rpplv equally to colieu nt 
type as well as non coherent type of search strategy The onlv difference between 
the two cases being the intei relation of parameters r d P r i uid P D for the ft p< of 
detector used This inter relationship is essential for computing the acquisition tune 
performance and should be deteimmed at the first instance 

3 7 1 Serial search acquisition parameters in the absence of 
code doppler 

In the modeling of serial search code acquisition piocess by signal flow gr iph method 
the cells, m code uncertainty region are represented as node md blanches between tlu 
nodes as transitions fiom one cell to anothei In such a flow giaph it is custom u\ 
[2 6, 16] to label each branch with the product of the transition probability assaulted 
with going from the node at the originating end of branch to the node at its terminating 
end, and an integer (including zero) power of a parameter denoted as » The parame tei 
a is used to mark the time as one proceeds through the giaph and its power repiesc nts 
the number of time units (dwell times) spent m traversing tint branch Fiutheimou 
sum of the branch probability (lettmg x = 1) emanating from each node equals unity 
The cells 3 are senally tested until it is determined that a particular cell corresponds to 

3 The term nodes and cells will be used interchangeably heie 
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the alignment of two PN codes within a fraction of a chip If tin number of code (hips 
to be seuched are N v and search proceeds in half chip lncrenu ut then the numbei cells 
m the code uncertainty region are q = 2 N u Only one out of these q cells represents 
synchronization l e it is assumed that conect phase lie in am one of q cells and sen il 
seaich can begin any where [2 16] The successive decision pioccss is presented In i 
signal flow giaph m Pig 3 5 A prion probability of achieaing coarse synchronization 
m the i fl cell pioyidecl it is not present in first (7 - 1) nodes is gnen as 


(9 -» + 1 ) 


(3 39) 


wlieie 9 — i + l is the numbei of remaining cells to be searched < ach possessing an equal 
probability of achieving synchronization Starting at cell 1 (Fig 3 0) if the synchronous 
state is detected the process will be completed m one dwell turn (r,/) as shown bv st itc 
r I11 7 tinnsfoim domain it is represented by Pd\Z weight of coi responding biiuch 
m Fig 3 5 If the synchronization has not been achieved the piocess will move to cell 
2 uicl tin corresponding branch weight will be (1 - Pdi)z Hue it encounteis titlici 
a false alum lcpiesenied by az K+1 ( a is false alarm probability and K is design 
paiametei to e\tend the integration time r d for verification of false alarm K < q) 01 
no false alum represented by (1 — ct)z At cell 2 the piouss will be repeated and 
search will pioceed fo check further cells until the synchronization is achieved In c ase 
synchronization is not obtained the search process wall start ill over again at stalling 
position 1 e cell 1 However now the overall probability of detection is increased in 
all the q cells toP m + P D2 {1 ~ Pm) and the whole search pioc edure will be repeated 
Since the SNR of synchronous state will be increased m each new retrace cycle aftei 
some passes synchronous state wall be detected and seaicli pioccss will stop 

Using flic st indaicl flow graph reduction technique [16] the transfer function of 
overall graph U{ ?) is determined This function is further use dm the derivation of 
a closed f 01 m analytical expression for mean acquisition time T acq and variance of 
acquisition lime a 2 acq The mean acquisition time is defined In [1 2 6 16] 

dlnUM ( 3 40) 

dz 



While the variance of acquisition time is defined as 



d 2 \nU(z) dAnUjz) 
dz 2 dz 


=1 


(3 41) 


Introducing tlie notation H(z) = (1 - a )z + az « « foi gam m advancing the se uch 
process fiom cell to (i + l)“ cell and ft as ove.all dotation probability m , ' 
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letiace cycle given as 


& = E P Dn n (1 - P Dm) 
n=l m= 1 


(3 42) 


Fioni signal flow gmph of Fig 3 5 loop gam of i th loop in j 11 k trace cade is gnen is 


(3 44) 


(3 43) 


L 3l = (l-p J )zH«- l (z) (3 13) 

Foiwaid pith gain of i ,h cell in j th retrace cycle is gnen b\ 

— p t (3jZ (3 44 1 

The tiansfei function of % th cell in j th pass is given as 

rp 

1-(1 -^)zH^(z) ( J ' 

Now liom Tig 3 5 the oveiall tiansfei function foi j 11 pass c in be written as 

Uj(?) = T 3 i + (1 - Pi)T^H{z) + (1 - P0(1 - P^TjiH^z) 

+(1 - P i) (1 - P 2 )(l - P 3 ) (1 - P q ^)T n H^\z) (3 46 

substituting Tp T) 2 , T 3 3 , , Tjq from (3 45 ) and Pi P 2 P q from (3 39) niU 

(3 46) and after some routine algebra the result can be expnsstd as 


(3 46) 


J [q 

which can be fuithei simplified i e , 


~~z~r r- A z) [- + -H(z) + -H 2 (z) + + -H^\z) 

- qqq q 


(3 47) 




(3 4b) 


Em N letiace cycles the oveiall tiansfer function of decision linking flov. graph isgi\i 




(3 49) 


Now using the definition of mean acquisition time from (3 40) and variance of atqm 
sit ion time fiom (3 41), the two parameters are found out m the appendices \ 1 uid 
A 2 respectively The final expression for mean acquisition time is 


-^L = ( q -l)(l + Ka) Ug-J-0oJ+l 


(3 3 0) 
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While foi vananee of acquisition time is 


a. 


acq 


T 2 

T d 


(q - 1)(1 + I<a) ^ 1 - ft, 


N 


N 

E 

j=i 


A 


3 + (9 — 2)(1 -T Koi) + 


Ac* (AT + 1) 
(g-l)(l + Ac*)J 


+ 


2{(g ~ i)(i + ft®)} 2 A fi -ft,\ 2 


a 


E 


ft 


(g- 1)(1 + Aa ) ^ 1 - /3j 


(g~l) 


(1 + Ace) 4- (1 + Ka) J 


N 


0 fq- j 


+ 


3=1 A? 
Ac*(A + 1) 


3 ill 


From (3 50) it is clear that mean acquisition time is a function of Pfa (denoted is a 
m the above equations), oveiall detection probability ft, number of cells in the (ode 
uncertainty legion q, integration time number of retiace tides N and false ilum 
penalty K If ftj is not too low only a few passes are required to iclneve synchiomzation 
The same can be said foi the variance of acquisition time However gnen the natuu 
and complexity of (3 51) it is difficult to draw meaningful conclusions regarding the 
mteiplai of vanous variables 


The mean and variance of the acquisition time for stand rid serial search appioach 
of Tig 2 5 me given as [2, 6] 


Tacq (2 — Pjd)(1 + Koi)(q - 1) 


acq 

Td 


2P d 


(3 j2) 


and 


vlcg _ | ( 1 + Ka ?q 2 (yz-p^ + pr) 

r d [ +6 q [K(Ii + l)a(2P D - P 2 D ){ 1 + Ka)( 4 - 2 P D - P 2 D )] + ^ 
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Hguie 3 5 Decision process signal flow graph for serial search code acquisi 
tion system of Fig 3 4 




Chapter 4 


Results and Discussions 


In tins chaptei we piesent the numencal lesults of the analysis earned out m the chaptei 
3 Firstly the results foi code phase detection performance of frved thiesliold method 
is piesented and discussed Secondly the case of adaptive thiesliold is considered and 
its results are presented The Nakagami m faded signal case is considered next and its 
detection performance is presented Finally, the corresponding results from the code 
acquisition point of view is piesented and discussed m the above order The manner in 
which these results are presented is suitable for making a comparative analysis with the 
existing known similar techniques and to clearly demarcate the resulting improvements 
and/or degradations 


4 1 Fixed threshold detection performance 

Tor a given P FA and N B the normalized threshold level f] can be obtained by nu 
meucal evaluation of (3 26) Upon substitution, this rf in (3 29) yields the P D foi 
that particular P FA provided the other parameters are suitably chosen The {P FA ) and 
(Pdi) arc fuithcr used to compute the mean and variance of acquisition tune m (3 oO) 

and (3 51) respectively 

4 11 Detection performance for unfaded signal in AWGN 

Fig 4 1 Using (3 26) P FA against the threshold level (rf) as a function of N B is 
plotted in Fig 4 1 With increasing N B larger threshold level n ' » needed tor a given 
value of P rA For a given N„ the rf increases with decreasing Pea 
Fig 4 2 The Pd, as computed from (3 29) vs SNR foi Pda = 01 Nb = 4 and 
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riyut 1 L Probability of false alarm Pp A vs normalized threshold level ij 

is i function of numbei of letiace cycle of recn dilation loop is plotted in Fig 4 2 In 
gene ill with hick lsiug ictiiee cycles the cmves shift tow uds the left and mcieise 
m detection piobibility is obtained Howevei m low SNR legion (< 0 dB) sufficient 
numbei of ictuce cycles are needed to achieve fairly high detection probability 
Fig 4 3 uiel Fig 4 4 The improvement factor in detection probability can llso be 
obt lined by men ismg N B) which is defined as ratio of detection probability for a given 
nurnbi r of post detection integration N B to the case when no post detection integration 
Np = 1 n employed The improvement factoi vs SNR as i function of numbei of 
leti icc cycles foi Pr i = 0 1, N B = 4 16 is plotted m Fig 4 3 A peak impioveme nt 
fact oi of 1 6 to 1 8 is obtained for N B = 4 and 2 5 to 3 1 is obtained for N B = 16 
Sunil u cmves foi Pj i = 0 01 aie shown in Fig 44 m which i peak improvement of 
2 9 to 3 2 time s foi N 13 = 4 and 8 5 to 9 2 times for N B = 1G is obtained As show n 
m Fig 4 3 and Fig 4 1 with decreasing false alarm probability highei impiovemeiit m 
detection piobibility is possible with increasing N B Moieoici the peak of the cunes 
shift towaids the left l e towards the low SNR region with lnei easing N B 


4 2 Adaptive threshold detection performance 

The threshold coefficient T for designated P FA foi gi\en shift legistei size M + 1 
and N b cm be obtained from numerical evaluation of (3 34) Upon substitution of 
T along with other parameters m (3 35), one can obtain the detection probability m 


m ... 
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Figure 4 2 Detection probability vs SNR for three retrace cycles 

eveiy letiace cycle of recirculation loop The mean acquisition time and variance ol 
acquisition tune then can be obtained by utilizing this Pfa uid Pd along with other 
parameters in (3 50) and (3 51 ) respectively 

4 2 1 Detection performance for unfaded signal m AWGN 

Fig 4 5 and Fig 4 6 The probability of false alarm Pfa against threshold coefficn lit 
T is plotted in Fig 4 5 and Fig 4 6 Tor specific N B and Pfa the threshold coefficient 
T dou cases with increasing M Moreover for wide range of Pfa bom 10" 6 to 10" 1 the 
variations in the threshold coefficient decreases with increasing M This is due to the 
fact that as more and more noise samples are included in the adaptive CFAR decision 
piocessor of Fig 3 3, the better will be the noise estimate and hence more accuiate will 
be the threshold coefficient, which in turn will results in accuiate threshold setting to 
achieve the desired false alarm probability 

Increasing the N e , . e , better estimation of decision variable Z, of Fig 3 4 the whole 
family of curves shift towards left as shown in Fig 4 5 and Fig 4 6 indicating again a 
reduction in T Moreover, for specific N B and T a reduction in P FA .s obtained with 

increasing M 


Chapter 4 


50 




riguie 4 3 Detection probability improvement factor vs SNR for four re 
trace cycles 

Fig 4 7, Fig 4 8 and Fig 4 9 The detection probabiht\ as obtained from the 
(3 35) for vanous SNR is shown m Fig 4 7 for designated Pr-i = 0 1, Nb = 4 and two 
diffeient M l e , M = 4 and M = 8 The similar curves wen obtained foi = 0 01 
and aie shown m Tig 4 8 Finally Fig 4 9 shows similar cuivc foi P rA = 0 01 N B — 6 
and M = 6 

The impoitaut observations are as follows 

• Due to iccuculation loop, the detection probability ciuves shift towards left 
ic, an improvement m detection probability is obtained m each retrace cicle 
The effect of recnculation loop is most visible m SNR lange —30 dB to +5 dB 
Beyond bdB all curves start saturating i e detection probability approaches to 
unity and below ~30dB another few retrace cycles are m c ded to obtain the fanly 
high detection probability 

• As shown in Fig 4 7, for given Pfa and N B further improvement is possible 
with increasing value of M in each retrace cycle of recnculation loop The same 
observation is also visible in Fig 4 8 

• As shown m Fig 4 9 increasing N B further improves the detection probability 









FALSE ALARM PROBABIUTY(Pfa) 
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« A comparison with Fig 4 2 and Fig 4 7 shows a vast impioiement m detection 
piob ability is obt lined by adaptive thresholding technique of Fig 3 3 

4 2 2 Detection performance for Nakagann-m faded signal m 

AWGN 

Fig 4 10 Tilt pdl of decision valuable when the lecened signal is seaeich f lelccl 
accoidmg to N ikigaiui m distnbution can be obtained bv (3 21) aie show n m Fig 4 10 
and Fig 1 11 foi vaiious fading paiametei rn The paianictu m < 1 indicates that 
u caved signal is moic severely faded than Rayleigh fading (m = 1) and m > 1 
com spend to tlu c isc s of less severe fading than than Ra\h lgh 

It c m be obsc ivc cl th at as m is mat vsed the symmetiv of c un es ma easts and the a 
be conic mon aid mon pc akc d Foi low SNR (—10 dB) peaking occuis at lngti aaluc 
of in while foi model itch high SNR (0 dB) the peaking occuis at compaiatiyely love i 
v due of m 

Fig 4 11 Tlu dc tcction probability as obtained fiom (3 38) is shown m Fig 4 11 
foi the pu mu tc is N a - i il/ = 4 Pr i = 0 001 and mac isiug \ dues of fidmg 
p name tei m As e an be obsoivcd in Fig 4 11 when the detection probability is m the 
i mge of 0 2 to 0 3, be (ta detection performance is possible with the it clue, mg fading 
seventy ( unit asmg w ) However below this iange of detection probability better 
de te t don pe i foi mane < with mci casing fading seventy (reducing m) The phenomenon 
is not unusual and has been observed by many authors [13 31] The reason behind 
t his phe noinc non is t h at the re exist an optimum N B for the fade d received signal below 
uul above which the performance degrades with mci easing m 

4 3 Acquisition time performance Fixed threshold 
method 

Fig 4 12 and Fig 4 13 The mean and variance of acquisition time foi modified 
serial s< auk slrategy of Fig 3 4 is computed from (3 50) Similarly the same par mi 
etas lot standard send search scheme of Fig 2 5 is obtamul from (3 52) and (3 i3) 
respectively Tin improvement factor for mean acquisition time is defined as ratio 
at mean auction tme when no Ration loop uproot to the mse roken tingle 
recnenMwn hop « present” A similar definition also applies to the improvemmt 

factor foi variance of acquisition time 


Chapter 4 


o5 


The mean acquisition time improvement factor for the case Nb = 14 16 number 
of cells in code uncertainty region q = 128 , false alarm overhead K — 64 false alum 
piobibility Pp 4 = 01 and 0 01 is plotted in Fig 4 12 and Fig 4 13 respectively 

An impiovement of 3 to 9 times is possible for mean acquisition time in the SNR 
langc of —40 dB to +10 dB Improvement factor for the variance is shown in Fig 4 13 
md 4 to 11 tunes impiovement is obtained With N B = 1 the peak of all the cunes 
Iks m the SNR region extending fiom -5 dB to +10dB while peak of the curves shift 
towuds left with mci easing Nb indicating that due to post detection integration maxi 
mum improvement is possible in extiemely low SNR region However rapid decrease m 
the impiovement factoi for mean acquisition time can be obseived beyond -5 dB with 
increasing Nb This rapid fall is due to the fact that with increasing N B recirculation 
loop increases the SNR quite rapidly while the detection piobabihty asymptotically 
mac iscs towards the unity md impiovement, finally settles down to a moderate value 
in be t vi c c ii 4 md 5 lh( same is true for variance improvement factor beyond 0 dB 


4 4 Acquisition time performance Adaptive thresh- 
old method 

Fig 4 14 The impi ovement factor for the mean acquisition time is plotted m Fig 4 14 
for the case N B = 1,4,8, q = 1024, K = 512 P FA = 0 1 and 0 01 The important 

observations aie as follows 

• An improvement factor of 4 to 13 is possible at P rA = 0 1 while 4 to 16 is possible 
feu Pi A as 0 01 in the SNR range of -40 dB to +10 dB 

. A c omp uison with the Fig 4 12 foi fixed threshold acquisition performance clear h 
shows better results are possible using the adaptne thresholding technique 

Fig 4 15 The improvement factor for variance of acquisition time is plotted 
Fig 4 15 for q - 1024, K = 512, M = 4, P F a = 0 1 and 0 01 An improvement factoi 
of 5 to 15 is possible, while similar curves for fixed threshold method m Fig 4 13 haae 
shown an improvement of 4 to 11 times 
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4 5 Acquisition time performance in Nakagami-m 
fading environment Adaptive threshold method 

Fig 4 16 and Fig 4 17 The improvement factois for tlu mean variance of the 
acquisition tune when (he received signal is severely faded is shown in Fig 4 16 and 
Fig 4 17 respectively loi vaiying degree of severity Foi fading more seeere (m < 1) 
thru R vylc lgh (m — 1) m improvement factor of 3 to 6 times in mean acquisition time 
and 5 to 6 limes in vanance has been obtained For Rayleigh fading case (m — 1) 
the two improvement factois are slightly higher In general as the fading severitv is 
reduced the improvement cuives shift upwards i e , higher improvement rs possible as 
the signal becomes moie and moie unfaded 



p(z*) 













Acquisition time variance improvement factor Acquisition tune variance improvement factor 
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Figure 4 16 Improvement m Mean of the acquisition time Nakagami-m 
faded received signal using adaptive threshold method 



F.gure 4 17 Improvement m variance of the acquisition time Nakagam.-m 
faded received signal using adaptive threshold method 



Chapter 5 


Conclusions and Suggestions for 
Further Work 


In this woik 1 icccntly proposed [15] approach for PN code icquisition based on the 
recirculation loop principle has been modified and analyzed m teims of code acquisition 
parametcis using fixed and adaptive threshold detection methods 

Pei foi mance ainlysis of these algorithms have shown that foi very low SNR lange 
(-AQdB to OdB) considerably high improvement is possible over the existing serial 
search techniques For model ate SNR range (OdB to lOdB) the improvement is still 
significant In particular, the average acquisition time and its variance can be reduced 
up to 5 timet, at SNR = -15 dB The overhead involved is not very significant as 
additional component e g post detection integration block (PDI) and adaptive decision 
processor cm be easily incorporated in the integrated circuits The foimei range of 
SNR is of mtucst dm mg the penod of initial synchronization 


5 1 Scope for further work 

Some of llie unie&olvtd issues axe 

• What should be the hltei bandwidth (B) ? ThB ans'wii depends upon the al 
located width ol frequency band of operation and paitially lies m the hands of 
frequency rcgulatoiy authorities and partially on the system variables such is 
the number of post detection integrations (. N B ) employed Since each test re 
quires summation over N B samples, the total search time will be proportional to 
N b P ) where P is the PN code period This implies that N B should be reduced 
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to decieasc the search time On the othei hand requnements of high detection 
probability and acceptable false alarm probability need Ng to be high This fact 
w is rlso obscived m chapter 4 in the case of detection pc iformance in Nakagami 
nr tiding environment Moreover, Ng = Br^ and B is chosen such that B — 1/T 
where J? rs simplmg peirod, therefore choosing faster sampling rate will increase 
B mcl hence Ng Now the choice of integration dwell time (t^) dictates the opti 
miz it ion pioblc m Thus parameters B Ng and tj. need c rrpful examination from 
opt inn/ it ion point of view so as to satisfy the constiamt of leduced search time 
on one hand and acceptable detection performance on the other For a given 
bandwidth B one possible approach as applicable to fixed threshold detection 
method could be is follows 

Foi i give n false alarm probability Bp a solve foi i] For a given de 
kction piolnbihty Pg l and SNR ( r y l ) solve transcendently for Ng now 
using flu ulation Ng = Bt^ determine the Ti Finally using the Td 
P s A md Pjji determine mean acquisition tune and its variance 

Sinnhi approach for adaptive detection method could ilso be devised keeping m 
view the fact that there is one more parameter M (number of reference cell in 
the ad iptive decision processor) that needs to be optimized larger value of M 
will give belt ci estimate of noise or interference Howewr, no procedure exist in 
the literature cited [18, 31] which gives optimum M foi given P FA and P Dl 

• Another fact or on which the acquisition time directly dept nds is false alarm penal 
ty (Iv) l c , number of integration needed to detect a filse preposition Though 
K should be chosen according to the criterion K -C (j procedure foi finding 
optimum Iv is f ir from being understood in contempoi uy literature 

• In the si uidaid serial search technique a single cock phase is examined at a 
time for possible synchronization A natural extension of this technique will be 
to use two or more branches foi simultaneous search of code phase uncertainty 
xogion This approach, on one hand, will reflect in the increase of hardware 
complexity, which is quite affordable with present state of the art ASICs On 
the oilier hand, by searching more than one code phase at a time it is expected 
that the acquisition time will decrease in direct relation to the number parallel 
branches used Specifically, devise a system where the entire q cell uncertainty 
legion is subdivided into N p > 2 equal components c ach containing ^ cells 
Each of N p branches will be responsible for searching one component m pai allel 
with remaining (N p - 1) branches Furthermore since the N p searches progiess 
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m paialld undci a common conti ol at any point along the way all N p alignrm lit 
ex imiiia) ions must be dismissed before search can proi ( < d further 

« When known code Dopplei is present the acquisition time performance of s\ s 
tom shown in hg 3 4 is iffected in two ways The first being lelative shift between 
loc illy gcnei itecl and leceived PN code phase to be tune varying i e relative 
phase between the two codes is smeared during the clwdl tune of mtegiate and 
dump Tlu it suit is eitliei an increase or decrease of detection probability de 
pending upon whether code Doppler is causing relative tode phase to mcreast or 
decieast Thus the signal flow graph of fig 3 5 should b< ippropnateh modified 
Tlu ‘second and potentially more dominant effect is that rode Doppler affects the 
avciagt seaich rate (sweep late), which m extreme oses can reduce it to zero 
tlu it by int leasing the acquisition time to a very high value and hence needs 
clitic al ox munition A bnef description for the effett of code Doppler on the 
ivti igt staicli latt is piovided m [2, 6] 

0 roi unknown tode Doppler, a code late uncertainty legion of the spreading code 
will be introduced which can also be partitioned into manv discrete bands foi 
cithci sanl oi pnallel search This search procedure will now be more compli 
cated m the sense that it is two dimensional now one in code rate and anothei 

m time 

0 The ideas ol parallel search and recuculation loop can also be applied to the othei 
loims ol senal code seaich piocesses such as sequential detection etc 
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Appendix A 


Derivation of Acquisition Time 
Parameters 


A 1 Mean acquisition time for serial search 


Taking iidtiu il loganthm of (3 49) we get 

<7-1 

In U(z) = - In q + In Y2 H l (z) + In 
i=o 


A 




AD 


differentiation of (A 1) with respect to z will give 


1 , - = gaogklhMk+lhf ±iLz Q]} 

dz ' J Et 1 
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dl 
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. At 
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dz 
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7=1 T j=i 
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dz 


[1 -(i-j8,)*£r*-»(*)p 

/I»“*(*)[(g - i){a(If + 1)A +1 + (1 - a)*} + (*)] (A. 3) 


Since JJ(z) — az ft *' 1 4 (1 - cv)z> theiefore, [if(z)]z=i = 1 so that 

d[zH q ~ 1 (z)\ 


[dl] 


dz 


= (g-l)(l + oA)+l 


(\4) 


z = 1 


Substituting (A 3) in (A 2) and recognizing the fact that S/Lo ^ iesu ^ 

aftei much simplification can be written as 


d[lnU(z)] 
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Since (q l)(i + An) >> 1 , therefore after some cumbersome manipulation the 
denied expn ssion for mean acquisition time is found out to be 

d[\nU( ■*•)]] 

= {q-i)(i + Ku) 


T, 


cuq 


Tl 


(b 


1 v- 1 1 


+ 1 


(4 6) 


A 2 Vai lance of acquisition time for serial search 


Stutmg with the double dilleientiation of (A 1), we have 


d 2 i tk \ d \ 
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Theiefou , van iuc ( of requisition time can be written as 
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(All) 


Substituting (A 8), (A 9) uid (A 10) m (A 11) and using the fret (g-l)(l + Ac*) » 1 
the icsult can be expressed as 


u acq 
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Appendix B 


Probability Density Function of 
Decision Variable 

B 1 Pdf of decision variable for no post detection 
mtegiation case Unfaded signal in AWGN 

Tian&fonmiig (1 ( )) (o the polu tooidmates le , w(i) = R(t)cob9 and v(t) — R(t)sm9 
To find tilt chslnlmUon of H and 0 w< simply transform (w,v) to (R 9) m two dimensions 
as the ic is oik to out (onespondtnot 

K v) — > (R, 9) 
p(u ), o)dwdv = p(R, 9)dRd9 
dwdv — RdRdd 



The random variable 0 is assumed to be uniformly distributed m the internal [0 2tt] 
theicfou distiibution of It can be found out by integrating tin (B 1) with respect to 6 
m the above mtcival 1 o , 
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WIkk /<>(') ls Bl( moclilucl Bessel luiulion of 6is{ kind ind /e 10 th oidei and is defined 
as[ll]i,(0 “ i C < M>( m>s0)r/0 

Smu ij(l) ~h (/) tlu k loti ipplymg llu huKhmuiiil fhcoiem of piobibiliti [12] cm 

(B 2) w< obi mi tin ( list libnl ion of ^(1) imdei signal present condition, (Hy 
potliesis U i ) l ( 


1 ( (—fir + 7*) I 0 fm y > 0 ^ 

ot ki wise 

nnl tlu 1 »U ol ^ (t) when the signal is absent, (Hypothesis H 0 ) can be obtained 
In substituting l — 0 m (B 3) mcl k cognizing the fict tint J o (0) = 1 





ioi y > 0 

0 C)thn\US( 


(B 4) 


B 2 Pell of decision variable for post detection in- 
tegration case Unfaded signal m AWGN 

LcUmg / I lu pdl ol 7 nuclei hypothesis Ho c in be wutlcn from (B 4) is 

, Na ( rY iN }i \ 

M>U) - \~~2^r ) d " 1 ( B °> 


To find mil I hi [nil ol dmsion vunbk Z , cliai ntenstic fuu< lion (CF) appionch [12] 
will hi followed u I ike I lie I ipluc ii insfoim ol po( y /) tdvc its Ng ih povei uul 
mvcit (lit Ii msloiiu to obt mi lh< pdl ol Z undoi both tlu lit po theses 
Ihcicloic uudci (lie hypothesis i/ 0 , Hit Cl denoted by il/o(s) is given as 


i\J 


0 




(~*N S \ 
V 2 ct 2 J 


exp(-s~,)f /7 


(B 6) 


(B7 > 

The c cju ilion (B 7) gives the CF ol decision vaaable Z unclu the hypothesis H 0 foi 110 
post elc U 1 1 ion ml c gi it ion c asc The icloie, the pdf of Z lor post detection mtegi ition 
case tan In c vpicsse cl as 


Po(Z) = Hj[M 0 (s)] N ‘«P( sZ )^ 


(B8) 
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^ /; (s + ^) % ( * y) 

Sm(( llui< IS isn.g. |).)1( of N/i lh oid( i U s = mdusuluc it this pole is found 
out <0 be 


Cn,~W zK ""^-^ 


(B 10) 


(lu k lot e in ikiiu, use of ( ' uu liy’s le sidue the oiem [14] the like lihood pdf of dee isioii 
M inblt / unda hypothesis II, is found end lolx 


/»..(/) { (^ / ) AJ ’ * (ir])7 (X P (~1?) ioiZ>Q 


(B 11) 


othcmise 


lo find out the pdf etf deeision vunble Z undei the hypothesis Hi \ie proceed is 
follows Sine < i Jf 1 " (lie ie foie (he pdf ol using (B 3) c in he expressed as 


e i N » 

P\\ i ) - } <M> 


iNh 1 r ( t 

M 2 v~ 


(B 12) 


I he (1 undet hypothesis lf\ denote el is Afi(s) is 1 hen de fun cl as 


4l/j (s )- j P\ U)exp(-W)^/ 


(B 13) 


uid e\p Hiding (lu modified Bessel fund ion [14] I 0 (i) = EEo^T^fr- AA(s) cm be 


expressed is 


AM-) 


^ N}j [* Njs^i ] ^ 


oo Ahn f 
r — v la 


(B 14) 


Re in uigmg (lie (e inis of y mel mlaehmgmg (he sign of nifegi ition and summation 


we hive 


;n(iNE 


r 4~rr- / ?i' 1 txp - TJ ^ 


(w») i '0 


(B Lo) 


The ml eye d J, e m be eviluated by making a substitution / (> + %*) - 9 and thm 
using (he definition of Gamma function [14] i c , T(n + 1) = «' = Jo° exp(-0) d $ 


(- + &)' 


(B 16) 



subsill ul mg lilt / 1 m (B i r )) m<l pciloimmg some loutuic algebiathe final lcsult <an 
b( expiessed is 




Mjl) 


- _ 

J 


7»* 

i' + t) 

- (xp 



There foie ptU of eletisum vmible Z untie 1 the hypothesis H\ is given 


as 


i> m) 


En 

i 


1 N„ 


2ni 


1 


(-+?) 




WN B 

s + S 1 


+ sZ 


d s 


multiplying m liiuud itoi mtt duioinmitoi of (B 18) by i\p (7, + J-)] 
kiting s | ^ - u so (hit f/s -- f/« tin p\{/) llun b( tomes 


(B 17) 


(B 18) 


■JV 

Ar || 

.2fT 2 . 

( \P 


Nn 7/ H 


Z 


1 J- 

(i 

( 1 

1 

JU 

+ 

t 

. 27r? y 


it 


du 


(B 19) 


Pin modilucl Btssd fuiulion of n lh oitlci c 111 be expressed 111 tiansfoim domain [11] 


r l 2 L 


I " ( " I ’ ) = (^J ^- J ^l^ 0 XV (l^ + T> ]d, 


(B 20) 


lolling a = \fl/ vud h ~ we e in wide 

l,N u l 


27T / 


i >„ 1 ')> ( / " 1 


1a‘Z\ N ‘~ l , IJZ%N 2 b 


u2cr 2 


du 


NhJ 


u -l I 2 [ 


2a 2 


(B 21) 


using (B 21) m (B 19) I lie likelihood pelt of lh{Z) cm Ikon be expressed is 

A n 1 

ft f r <s "1 **■ 

■ « - (r 


pit 


exp 


A«(7, I |i) 


2Za 2 


HNb j 


i«, -1 I 2\/ 1 (B 22 ) 


2er° 


Fin illy all e 1 some routine ilge In 1 flic the result can be expie ssecl m simpler form 
Pi(Z) = 


Nn-I 




0 


-Nn ( 7 * + ^x)] [ 2 ^ 7 ] 2 j ^b-i ( 2 V^Fj fol Z - ()b 23) 

otherwise 


Normalising Z by 2a 2 /N B = N 0 r d or, equivalently, letting Z = we can reiuite 
(B 11) and (B 23), re spee lively, m 1 simpler form 
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uicl 


P (^) = 


{/ ) N n-i , „ x 

W-L)>“^P(-^ ) 
0 


foi Z -> () 
ot heivus( 


(B 24) 


M(^) = 


V W;i7 
0 


<\p(-Z H - iVflT ,)In b -i{2\/N]P) Z ) foi Z >0 

otheiwi&e 


(B 27) 


B 3 Pdf of decision variable for post detection m- 
tegration case Nakagami-m faded signal in 
AWGN 


1 xpussmg (lit luodilu (1 Bt ss< 1 lime lion of (N B - l) tk oidn m suies sumnntion ionn 

IS [1 1) 




pulling (B 2(>) m (B 26) Hit k suiting funclion will be 

EilJ. 

00 

ud-z'-i.jvdE 




(N„Z 


,=0 > -l ) 1 

Substituting (1 1 ( )) uni (B 27) m (3 20), lli< iesulting funclion will be 




n ; 


/ 


lilt -t 




■d7» 


(B 26) 


(B B 


(B 28) 


0,,y n . 7i "ozi oo 7 

,,||/ ' 1 iwfe) "^-^Ster ,-ijf 


00 


7 , m ~ 1+7 exp 


f m , T 
-7» I — + v /j 


n 


d h 


h 


(B 29) 


lolling 7, | jV/d - 0 <uul using tlie definition of T( ) function [14] the integial I\ 

it duces to 


/, = 7 — T7q7 r e m+r - l cxp(-B)de 

(# + "*) /0 


r(ni + r) 


(B 30) 


\ppcudi\ B 


Subsiil u( (13 30) m (13 29) (lit i ( suit will be 


, /M (m/il)'" \ 7 ' 00 (N 7 ,iu±_zl 

Pi{/ ) ~JT m ) A r f ^_L I> + 0 

( ’ lA, " J 'a,T (W b+.j^+w,)* 


(B 11) 




nr f f/ 1 ] i ' rrl np(_7MV'^^ r(w + l) 

1 <"'>(«' "»)' ,( Z) S 75 vM 0 ^f (B,2) 


k tlmg 11 hI” Ulnt ll ’B ( mi}-, — i^(7i) uid leiu uiging Hu t-cmis \ie get 


{^) Nli 1 ( tp(-Z ) T(m + j) f uZ 
l(m) (l + !/)»» “ i T(N b + ? ) {T+h 


(B 33) 


Mulliph mg m umiK t Hot uid diuoumi itoi In V(N B ) tin find u suit on b< cxpussul 


lh(/) W* 'until T ( M uZ 
(i i u) m [N }i - i)i lJl \ Nn T+7i 


(B 11) 


' ' ' " / ' - i / 

VVIhk ,/,(« I) i) is confluent Hypei geometric function defined is [3 14] 


, , , V 1 ( 6 ) I*r(# +a) 

' y|( '''’' ,) = rw?TWTA 


(B 3j) 



Appendix C 


Evaluation of Detection 
Probabilities P FA and P D 


C 1 Unladed signal m AWGN Fixed threshold 
case 


C 1 1 Pi obability of false alarm for post detection case 


rill piob dulil's <>l hist il uni Pi \ for fixed signal in AWGN can be obtained In 
mligiiling Ilu likelihood pdl P 0 (/ ) 



P, 


i i 


Po(Z*)dZ 


{7*) Nl > l cxp (-£*) 
(N„- J)' 


t IZ 


Using lh< ddunlumol nuotnpldi g unma fund ion k , 


I {n /) r 0" 1 c\pH)t/(9=(7i-l)'(xp(-o ii4 
I « j>= o P 


(C 1) 


(C 2) 


Using (( 2) ui (( I) Uu ksuK i m In (\pnssccl is 


Nrt-l „+p 

Pi a = cxp(-/f) £ ~r 

p~Q P 


wliuc //* is nmmaliAcl tlncshold and is defined as rf ^ 


(C 3) 
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c 1 2 Pi obabilit y of Use alai m foi no post detection integi a- 

(ion case 

il„ piob lbililA ol i >K( il inn uudei no post detection integration case c in l, f 
obi nut (l In subsl tfulnu, Nn I in < qinhon (C 3 ) i c 

Pi i-=(\p(-//') (C 4) 


C 1 3 Pi ob ability of detection for post detection integration 
c asc 


flu piolnlnlih of d< (( t lion P, h lot ulna cydt ol it( nail ition loop can be oh 
I unt (1 In 

Pd,= r M Z )( 1 Z* 

>‘l 

l,U /, (v/oj 1 mK-/'* - )dZ* 

1 

fyM <M>( 9 , N n ) (C i) 

rill liiodilu d Ibssd 1 1 tin I n >n i( ) t m lx rypitssed in stilts sum foim i c 




Ia±_ 


,=o ''(«•„+, -if (t6) 

(Substituting (C (>) m (( 9) uul ptiloimniR stunt loulmt ilgduatlie ltsuit cm lx 
(\pnssttl is 


Pn, <\p( vV/*)£ 


(Nn^Y 


rgiWH » -i)'& 


00 


(Z') 




:\p {-Z*)dZ* 


(C 7) 


/V) 


Ilu mli gill /v i tan ix t\ dual cd using llic dc fimtion ol mt omplete g imnn func l ion 
given in (C 2) i t 


Nb+i-I ( n )l 

A i - (A^/H t - 1)'< cp{-7]*) E ”T (C8) 

p = 0 

finally subsl if utmg (( 8) m (C 7) the d( ted ion probability expiession can be expressed 


as 
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SO 


P D = exp(-T N„-n) £ ' °t ' 1,1 y 


r=() 


/=() 


P' 


(( 0) 


C 1 4 Probability of detection for no post detection integra- 
tion 


The detection piobibihtv foi no post detection liitegiation ( isc ma\ be obt uned In 
putting N b = 1 in equition (C 9) l e 

00 V In V 

Pdl = exp(- 7l -77 i() ) 

=0 1 p= 0 p 

Which if desned cm be c\piessed in Marcum Q function hum 


C 2 Unfaded signal in AWGN Adaptive threshold 
case 

C 2 1 Probability of false alarm for post detection integration 
case 

The piob ability of false alarm foi a given -value of T is obt untd by mtegi iting the 
detection piobabilitv conditioned on a tlneshold rj =T\ mu p(\) 


f 00 f oo 

Pp± = p (\) / Po (Z*)dZd\ (C 11) 

Jo Jt\ 

Since \ is defined is 

JU 

\ = 52 Z, foi (C 12) 

t=-J£ 

wheie each of Z, consists of N B independent white Gaussian noise samples the pdf of 
X therefoie cm be written as 


p(V = 


\ L 1 exp(— A) 


where L = M V B 


(L~ 1)> 

Thciefoie substituting (C 13) m (C 11) we have 

\ L ~ l exp(- \) f°° (Z ) jv b~ 1 ) e\p (-Z ) 


Ppa 


(L-iy 


(Nb - l ) 1 


dZ d\ 


p (z ) 

" 1 V 

h 


(C 13) 


T\ 


(C 14) 
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SI 


using the definition of incomplete gamma function gnen in [1 i] the mini lntt^i »1 / j 
can be eaaluitc d 1 e 


/i = (Nb — l) 1 trp(—T\) 


N b - 1 

E 


7=0 


(T\y 

p' 


(C la) 


substituting (C 15) m (C 14) and changing the oidei of summ it ion and mte 0 iatioii \\c 
get 

i 

Now mal mg use of definition of Gamma function [14] the mtcgi il I> is found out to 
be (L+p — l) 1 Substituting it in (C 16) and peifoimnig sonn nnnoi algebu tin fin il 
icsult can be wntten m simplified form i e 


P 


F 1 


1 n ^(MN b +P- l\ T> 
(1 + T)U»b ^ 


V 


IB 1 -p 

V 


/ 


(1 + T)P 


(C 17) 


wlieie notation ^ J is used to denote when a and b uc positm mtcgei uul 

also foi repiesentmg ^ hen 1 ind b are P osltn( 11011 mte S (r 


C 2 2 Probability of false alarm for no post detection mtegi a- 
tion case 


the piobabilitv of false alann foi no post detection mtegiation case ma\ be obtinu d 
ba substituting Nb = 1 in (C 17) 1 c 


Pfa (1 + T)* 1 


(C IS) 


C 2 3 Probability of detection for post detection integration 
case 

The detection probability foi a given value of T is obtained In averaging the detection 
pxobabiht} ovei the pdf p(\) le , 

A L ~ 1 exp(-\) [«■ / , xphz . _ %Nb)Ini _ i( 2 ^N b Z 7 ,)dZ d\ 

J 0 (L-l)' JTX __ , 

Pi(z ) 


l 


(CIS) 
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Expanding the modified Bessel function /yv B -i( ) into sen< s suiniii ition fonn 1 c 


n.-i(2 Jn b z ->) = £; 

r=0 


(iVaZ Q )^— 

i'(N b + i - 1)' 


(C 20) 


Substituting (C 20) into (C 19) and intei changing the sign of summation and mtegi i 
tion we lii\t 


00 / AT \? 

S ^ + >-1)1 eM3( " '' N ° \P Z )dz (C >1 ) 

l 

The mtegial L, can easily be evaluited using the definition of incomplete gamma time 
tion [14] i e 


h = T(N b + i T\) = (N b + 7 - iyeip{~T\) V (C 22) 

,zr. p' 

putting (C 22) into (C 21) and aftei some simplification the i< suit can be written is 
_EEL (iY R a,V Nb+i-i j >p r oc 

p o» = £ r /r _ iV ^p(~7 ? Vb) E Tf/ \ I "' p - 1 exi,{-(l + r)\}ei\ (C 23) 

r=0 1 \ J - 1 x ) „=n P y° 

' — ^ 

/ 1 

The mtegial 7j can be evaluated using the definition of gunnn function le 1\ — 
(L+p—iy upon substitution m (C 23) Finally aftei sonic e umbeisome mampul itiou 
the lesult can be expiessed is 

D _ exp (~^N b ) °° (N Blt y Nb +A~ 1 ( MN b +P-l\( T \ p , r , M 

Pdi ~ (i TTF^ r ? 0 -7r- g ( p JlTErJ (C24) 


C 2 4 Probability of detection for no post detection integra- 
tion case 


The detection probability foi no post detection mtegi ition e ise mat be obtained In 
simply putting Nb = 1 m c quation (C 24) i e 


C 3 


P _ cxp(- 7l ) A ( M+p-i \ / T \ p 

Dl (1 + Ty<tor'h\ P JU+T) 

Nakagami-m faded signal m AWGN 
threshold case 


(C 2b) 

Adaptive 


to o roo 

P Dl = / p(X) Pl (Z )dZ d\ 

Jo Jt\ 


(C 26) 
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(°° \ I 00 Z ^ B ~ 1) exp(-Z ) _ / . uZ \ 1 „ 

' p(\) / 7- r — 7-77 771F1 771 A a ) r/Z d\ 

0 iV ; [h \ (1 + u) m {N B - 1) 1 V l + it) 


(C 27) 


m + r-\\(uZ V ,vl „ ,r ,0 

/, p(x) L lT T uym- 1)' S , ) IttH f,z (f S) 


1 ^ / rn + r - 1 \ (il;) 


T _ x V Vj±HZ / (7 \AjH- 

(l + «) l r ^o\ I J {Nb + 1 -D'I t^ ... 

I 

The mtcgn.1 I 2 is incomplete gamma function and is defined is 


/ OO 

/ (Z ) A » J - ~ l fip(-z )dZ (C 29 ) 
/r\ 


r oo * 

r(n 9) = / e~H n ~ l dt = (n - l)if® V 

^0 > — n 


(C jO) 


Substituting C 30 m C 29 we get 


e~ TX ~ / m + r - 1 W it V Aj ^" 1 (T\)* 

(1 + tt)™ \ 7 )\l+u) fa, p' 


P -T\ 00 

/j = 7 — — 

(1 + ») m r =0 


K 11) 


Now p{ \ ) is given b\ 


p(M = 


\ L 1 eip(— \) 


^ ' (L — 1)' 

substituting (G 31) and (C 32) in (C 29) we get 

,00 \i-lp,p(-\) e _TX A / m + T-l\ / <7 y^+r-l ( r x)»> 

L ! T 1 M 7 I 1 .Am ^ l _ ) l 2 4- (/ / “ nt 


(C 12) 


(L-l)' (1 + «)*»£& 


X “|m + r-l| ( u \ Sn fr l T*_ 

(1 + u) m (L — 1)* y r Ju + J fa p' 

f°° \ L+P ~ l evp(-(1 + T) \)d\ (C 34 ) 

Jo 


(C 13) 


substituting (1 + T)\ = 9 the integral J 3 reduces to guiiuii function defined as 
T (n) = / 0 °° 9 n ~ l e~ e d6 Hence integral 1% reduces to 

/ = r (F + P)_ (C 33) 

3 (1 + T) L+p 

Therefore Pdi could be written as 


1 “/m+r-lN/ u y"^- 1 T>(L+p-l)' ,, ,,, 

P °‘ ~ (1 +u)“ 2 i r ju+ttj p'U - 1)'(1 + r ) i+p 



\ppench\ C 


S4 


Fmallv leaningmg tin like poyyei terms we get the desired expiession fox detection 
probability 1 e 


Pdi = 


1 


(i + //)"(i + T) 1 1; 


E 


m + r — 1 


1 ii 


E 

i=o 


{ i+p-n 

( 1 \ 


U + 7 ) 


(c .7) 


C 4 Rayleigh faded signal m AWGN Adaptive 
threshold case 


The probability of detection foi Rayleigh Faded signal in C V can be found oul m 
n mannoi sinnln to Nnkagaim-m faded case Since Rayleigh f icling is the spin d c xsc 
of Nalugami m fading theiefoie dnect substitution of m = 1 yyili yield the detection 
pi ob ability expiession i e 


Pd = 


i ^ f_u V Nb ^T 1 f L+p-l 

((1 + u)(l + T) x ,_ 0 VI + m/ p=0 1 p 


) (iTt)' (c 
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